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Abstract 

We compute the cosmic relic (dark matter) density of the hghtest supersymmet- 
ric particle (LSP) in the framework of minimal N = 1 Supergravity models with 
radiative breaking of the electroweak gauge symmetry. To this end, we re-calculate 
the cross sections for all possible annihilation processes for a general, mixed neu- 
tralino state with arbitrary mass. Our analysis includes effects of all Yukawa cou- 
plings of third generation fermions, and allows for a fairly general set of soft SUSY 
breaking parameters at the Planck scale. We find that a cosmologically interesting 
relic density emerges naturally over wide regions of parameter space. However, the 
requirement that relic neutralinos do not overdose the universe does not lead to 
upper bounds on SUSY breaking parameters that are strictly valid for all combina- 
tions of parameters and of interest for existing or planned collider experiments; in 
particular, gluino and squark masses in excess of 5 TeV cannot strictly be excluded. 
On the other hand, in the "generic" case of a gaugino-like neutralino whose annihi- 
lation cross sections are not "accidentally" enhanced by a nearby Higgs or Z pole, 
all sparticles should lie within the reach of the proposed pp and e+e" supercolliders. 
We also find that requiring the LSP to provide all dark matter predicted by infla- 
tionary models imposes a strict lower bound of 40 GeV on the common scalar mass 
m at the Planck scale, while the lightest sleptons would have to be heavier than 100 
GeV. Fortunately, a large relic neutralino density does not exclude the possibility 
that charginos, neutralinos, gluinos and squarks are all within the reach of LEP200 
and the Tevatron. 
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1. Introduction 



It is by now well established that the observed, luminous matter in the universe cannot 
account for its total mass. Cosmological mass densities are usually expressed as ratio 
Q = p/pc, where p^ — 2 ■ lO^^^/i^g/cm^ is the "critical" mass density that yields a flat 
universe, as favoured by inflationary cosmology 0; p < (>) Pc corresponds to an open 
(closed) universe, i.e. a metric with negative (positive) curvature. The dimensionless 
parameter h is proportional to the Hubble "constant" H describing the expansion of 
the universe: H = lOOh km/sec Mpc. Observations yield 0.5 < h < 1. Even if one 
broadly defines luminous matter as everything that emits any kind of electromagnetic 
radiation, when averaged over the volume of the (visible) universe it cannot give Q > 0.01. 
In contrast, from the observed orbits of hydrogen clouds around a variety of galaxies, 
including our own, one derives |l| > 0.1; and from the motion of (clusters of) galaxies 
within superclusters one can deduce [||] f2 > 0.3 — 0.4. Finally, as indicated above, $7 = 1 
is predicted |2[ by models of inflationary cosmology; such models are currently favoured, 
since they can solve other cosmological problems, e.g. the flatness, horizon and magnetic 
monopole problems. 

The nature of the missing or dark matter (DM) cannot be derived directly from present 
observations. However, within standard (Big Bang) cosmology, the observed abundances 
of light elements (D, He, Li) can only be understood if the total baryonic mass density is 
less than about 0.1 pc- The discrepancy between the lower bound on the total mass density 
and the upper bound on the baryonic one has led to speculations that some neutral, 
weakly interacting stable particle might provide the bulk of the mass density today. In 
particular, it has been observed almost 10 years ago that the lightest supersymmetric 
(SUSY) particle (LSP) is a good candidate for dark matter. Its stability is guaranteed by 
a discrete symmetry called R parity, which can be imposed in most phenomenologically 
viable SUSY models; in particular, R parity is automatically conserved in the simplest 
realistic SUSY model |^, the minimal supersymmetric standard model or MSSM. This 
explanation is especially attractive since the primary motivation for SUSY has nothing 
to do with the DM problem; rather, SUSY automatically provides a DM candidate "for 
free". Moreover, for "natural" choices of parameters (to be specified below), Q turns out 
to be p of approximately the right order of magnitude. 

The primary motivation for the introduction of supersymmetry stems from the obser- 
vation [§] that in SUSY models, large hierarchies between mass scales are automatically 
protected against (quadratically divergent) radiative corrections, in contrast to non- 
supersymmetric models. In particular, within the non-supersymmetric standard model 
(SM) it is extremely unnatural to assume the scale of electroweak symmetry breaking, 
characterized by the vacuum expectation value (vev) of the Higgs field (H) = 175 GeV, to 
be much smaller than the scale of Grand Unification Mx — 10^^ — 10^^ GeV or the Planck 
scale Mp ^ 10^^ GeV. On the other hand, in the MSSM radiative corrections are under 
control, provided only that the mass scale of the superpartners is not much bigger than 
1 TeV. The relation between the scale of electroweak symmetry breaking and sparticle 
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masses is even more direct in minimal supergravity (SUGRA) models where the break- 
down of electroweak gauge symmetry is induced by (logarithmic) radiative corrections to 
the parameters of the Higgs potential [llO| , |^ . These models also have the practical advan- 
tage that they allow to describe the whole spectrum of superparticles (sparticles) in terms 
of a small number of free parameters. Moreover, the resulting sparticle spectra almost 
automatically satisfy constraints from K and B physics [11|, and lead to small additional 



contributions to electroweak observables [|T2| , in agreement with LEP results [jl3| . Finally, 
precision measurements at LEP have shown [|l4l that the non-supersymmetric SM does 
not lead to a Grand Unification of the gauge couplings, whereas in the MSSM all three 
gauge couplings meet at scale Mx — 10^^ GeV. While this result does not depend |]15|, [16| 
on the constraints on the sparticle spectrum imposed by minimal SUGRA, it does lend 
credence to the assumption that there is no additional threshold between the weak scale 
and Mx] this assumption is an important ingredient of SUGRA models with radiative 
symmetry breaking. 

Quite a few papers have already been published in the last decade that contain cal- 
culations of the relic LSP density in some version of the MSSM. However, older papers 
often assume rather light sparticles, in conflict with recent experimental bounds. 
Moreover, many previous calculations ||^, l7, 18, 19, 20 1 involved simplifying assump- 



tions about the sparticle spectrum, which were often in conflict with SUGRA predictions 
(e.g. by assuming the LSP to have the same mass as squarks). Some computations 



2J] do take the SUGRA relations between sfermion and gaugino masses into 



account, but the additional constraints on model parameters imposed by radiative gauge 



symmetry breaking have still been ignored in these papers. Moreover, in refs.|21, 22, 24] 
only the case of a "light" LSP, with mass below that of the W boson, has been treated. 
Indeed, only two calculations of the annihilations cross sections of a heavy LSP exist to 
date llTSi, [T^, and neither of them is fully complete. In ref. [O] the annihilation into one 



Higgs boson and one gauge boson has not been included; moreover, the given expressions 
do not seem to be applicable if the mass of the LSP is less than half the mass of the 
heaviest Higgs boson. Ref. |19[ does treat the full list of possible final states, but only 
considers unmixed (pure) neutralino states; we will see that this actually gives wrong 
results for the gaugino-like LSP even in the limit of infinite LSP mass. The results of 
ref.[ff9ll have been used in refs.[EOl |23|, |25ll as well. 



|25|1 as well. We have independently computed all 
relevant annihilation cross sections for a general, mixed LSP eigenstate. This part of our 
work should be useful beyond the context of the more restrictive SUGRA models. Very 
recently another calculation of the relic density of heavy LSPs has appeared The 
authors state that they generalized the results of ref.|T9| by including neutralino mixing, 
but no explicit expressions are given; moreover, no SUGRA mass relations are assumed. 

We are aware of only two calculations of relic LSP densities ^ in which the 
constraints imposed by radiative gauge symmetry breaking have been taken into account. 
However, in these papers a specific SUSY breaking scheme, the "no-scale" ansatz PH| , 
has been assumed; in this scheme the LSP density turned out to be too low to account 
for all dark matter.S Moreover, the Higgs sector has only been treated in the tree-level 
approximation. The importance of Higgs exchange contributions has been pointed out 



^^This result can also be derived from refs.[^2ljB3] once one makes use of the fact that these models 
cannot support a higgsino-like LSP; see also ref. ||24| . 
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in ref.[^, again using tree-level formulae for Higgs boson masses. However, radiative 
corrections to these masses can be [^] substantial if the top quark is heavy, rrit > 100 
GeV, which now seems likely.! 

In this paper we compute the LSP relic density in minimal SUGRA models with 
radiative gauge symmetry breaking. We use 3 free parameters to describe SUSY breaking, 
as opposed to only 1 in refs.[^ Due to the constraints imposed by radiative gauge 



symmetry breaking, which in particular fix the Higgs spectrum for a given set of SUSY 
breaking parameters and given m^, the resulting parameter space is still sufficiently small 
to allow for an exhaustive scan; we thus do not have to rely on "simplifying assumptions" 
of often dubious validity. As indicated above, we always include neutralino mixing, as 
well as all kinematically accessible final states, when computing the LSP annihilation 
cross section. Our analysis of radiative symmetry breaking includes effects of the Yukawa 
couplings of the b quark and r lepton, which can be quite important Mixing 
between the superpartners of left and right handed fermions is also treated exactly. One- 
loop corrections to the Higgs sector are included, and all experimental constraints on 
sparticle masses are taken into account. 

We find that the model can easily yield sufficient dark matter to close the universe, if 
the SUSY breaking common scalar mass m at scale Mx exceeds 40 GeV. This conclusion 



is closely related to the result of refs. ||33| , [3^ that a single light slepton suffices to make the 



relic density of a gaugino-like LSP uninterestingly small. A qualitatively similar result has 



been found in ref. p4[ in a more general context. On the other hand, requiring the LSP not 
to overdose the universe does not lead to upper bounds on sparticle masses which are both 
relevant for existing or planned experiments and are valid for the entire parameter space. 



Of course, it is quite unnatural |BB|, ITO] to assume sparticles to be much heavier than 



the W and Z bosons, but naturalness arguments cannot be translated into strict upper 



bounds; even a rather high [T^, [T^ bound from cosmology would therefore have been 
welcome. Alas, there are two different ways in which such a bound can be circumvented. 
One possibility is to have a light higgsino-like LSP, with all SUSY breaking masses being 
very large. Note that the higgsino mass does not break supersymmetry; moreover, we 
argue that in such a scenario, supersymmetry would be extremely difficult to discover in 
laboratory experiments. The constraints from radiative gauge symmetry breaking exclude 
the possibility to have a higgsino-like LSP if the top quark is heavy, the precise bound on 
rrit depending on the ratios of the soft SUSY breaking parameters. The second possibihty 
to allow for a very heavy sparticle spectrum is to choose the LSP mass to be close to half 
the mass of the pseudoscalar Higgs boson; this strongly enhances the annihilation of the 
LSP into SM fermions, in particular b quarks and r leptons. Within the framework of 
minimal SUGRA, this scenario necessitates a large ratio tan/5 of the vevs of the two Higgs 
fields of the model, but such solutions can be realized quite easily |31|, |3^ . 

The remainder of this paper is organized as follows. In sec. 2 we briefly describe the 
formalism necessary to compute the DM density for a given LSP annihilation cross section. 
We also list all possible LSP annihilation processes, and give a short description of the 
parameter space of minimal SUGRA models. In sec. 3 we present illustrative examples 
of LSP densities in such models. The effects of s-channel poles as well as thresholds. 



^Leading one-loop corrections to the Higgs sector are included in ref.|gj|, but only for light LSPs. and 
without the constraints imposed by radiative gauge symmetry breaking. 
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where new annihilation channels open up, are discussed. We also explicitly demonstrate 
the importance of the SUGRA imposed running of sfermion and Higgs boson masses. 
Finally, we study the dependence of the DM density on the free parameters of the model 
by means of several contour plots. Sec. 4 is devoted to a discussion of the bounds that 
can be derived from the requirement that the LSP relic density lies in the cosmologically 
interesting region. As already mentioned above, a strict upper bound on sparticle masses 
can only be derived if we artificially restrict the parameter space of the model. In contrast, 
a non-trivial lower bound on slepton masses can be derived from the requirement that 
the LSP relic density be close to the critical density; however, the masses of the gluino, 
the light chargino, and the squarks can all be near their present lower bounds. Finally, 
in sec. 5 we summarize our results and draw some conclusions. Complete lists of all 
LSP annihilation matrix elements are given in Appendix A, and Appendix B contains an 
example of a check of some amphtudes for the production of longitudinal gauge bosons 
using the equivalence theorem. 

2. Formalism 

In this section we describe the formalism necessary to describe the numerical results of 
sees. 3 and 4. We first (sec. 2a) give a short summary of the calculation of the present day 
DM density for given mass and annihilation cross section of the dark matter candidate x- 
Sec. 2b contains a brief description of the MSSM; it also contains a list of all annihilation 
channels of our DM candidate. The resulting annihilation cross section will turn out to 
depend on the whole sparticle and Higgs boson spectrum. In sec. 2c we therefore give a 
brief summary of minimal SUGRA models, where the whole spectrum can be computed 
in terms of 4 free parameters. 

2a. Calculation of the DM density 

We begin with a brief description of the calculation of the present relic mass density of 
a DM candidate assuming that the mass as well as the annihilation cross section 
c"ann(xX anything) are known. Following the prescription of refs.0, |T^, we first intro- 
duce the freeze-out temperature Tp, below which the XX reaction rate is (much) smaller 
than the expansion rate of the universe. It is convenient to express (inverse) temperatures 
in terms of the dimensionless quantity x = m^/T; the freeze-out temperature can then 
iteratively be computed from 

, 0.07Q4Mp(a + 6b/xF)c(2 + c)m^ 
XF = m -. (1) 

y/9*XF 

Here, Mp = 1.22 ■ 10^^ GeV is the Planck mass, and a and b are the first two coefficients in 
the Taylor expansion of the annihilation cross section with respect to the relative velocity 
V of the XX Pci'ir in its center-of-mass frame: 

V ■ (TannixX anything) = a + bv"^; (2) 

notice that v is twice the velocity of x in the xX cms frame. Furthermore, g^, is the effective 
number of relativistic degrees of freedom at T = Tp. A highly relativistic boson (fermion) 
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contributes 1 (7/8) to g^, whereas very nonrelativistic (slow) particles do not contribute at 
all. However, often Tp turns out to be close to the mass of the heavy particles of the SM 
(the c, b or t quark, the r lepton, or the W or Z boson); in this case a careful treatment 
of the threshold is necessary |^ if ugly jumps in the curves are to be avoided. Finally, the 
constant c in eq. (|lD is a numerical parameter introduced IQ to achieve smooth matching of 
approximate solutions of the Boltzmann equation above and below Tp; following ref. 
we chose c = 1/2. Given xp and cTann, one can compute 



2_ Px _ 1.07- lOVGeVxi. 



''""^ - pjh^ V^Mp (a + 36/xi.) ' ^ ^ 

where the rescaled Hubble constant h and the critical density Pc have already been intro- 
duced in sec. l.B 

Eqs.(P-(0) describe a simple, approximate solution of the Boltzmann equation that 
determines the abundance of any particle species. While not strictly correct in most 



cases this treatment reproduces the exact numerical solution to 10-20 % accuracy || 
given that h'^ in eq.(|^) is only known to a factor of 2, this accuracy is fully sufficient for 



our purposes. However, it has recently been pointed out that there are three cases 
in which this approximation fails badly: close to a threshold where a new annihilation 
channel opens up that dominates the total annihilation cross section; close to a very 
narrow s— channel resonance; and if the next-to-lightest sparticle x' is close in mass to 
the LSP, and dannixx' ~^ anything) ^ CTannixX ~^ anything). Unfortunately, in these 



three cases the proper treatment ||3^ is considerable more cumbersome than eqs.(||)-(|^). 



For reasons of computational simplicity we will therefore use the approximate treatment 
throughout, but we will be careful to point out the situations where it might fail, and will 
qualitatively describe the result of the proper treatment in such cases. 

Note that xp in eq.(|]) is almost independent of and (Jann] we find 15 < xp < 30 
for experimentally allowed choices of parameters (see below). The number of degrees of 
freedom g^, at Tp ^ m^/20 increases monotonically with m-^, but again the dependence 
is quite small: 8 < ^/g^ < 10 for 20 GeV < < 1 TeV. The details of the model (particle 
masses and couplings) affect the prediction for therefore predominantly through the 
coefficients a and b describing the annihilation cross section; we now briefly describe the 
calculation of this cross section. 



2b. The annihilation cross section 



Within the MSSM, only the lightest neutralino is left |]38[ as a viable DM candidate. A 



stable LSP has to be |^ both electrically and color neutral, since otherwise it would have 



been found in searches for exotic isotopes. This leaves one with the lightest neutralino 
and the lightest sneutrino, but the latter is excluded by a combination of the bound on 



* Within the framework of inflationary cosmology, eq.(y) should strictly speaking be interpreted as a 
calculation of the expansion rate, rather than a calculation of the density. After all, inflation predicts 
j2j f2 = 1 to high precision, independent of the details of the particle physics model. However, only for 
a small range of values of the (absolute) mass density p does H, and thus the age of the universe, come 
out close to the observed value. 
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the invisible decay width of the Z boson obtained at LEP [40| and the hmits derived from 



the unsuccessful search for relic sneutrinos using Germanium detectors ||41|| . 

The sparticle spectrum of the MSSM contains [§,0] 4 neutralino states: The superpart- 
ners of the B and VF3 gauge bosons, and the superpartners of the neutral Higgs bosons 
li\ and H\ with hypercharge -1/2 and +1/2, respectively. However, after electroweak 
gauge symmetry breaking these current eigenstates mix; their mass matrix in the basis 
{B,W^;h\;hf) is given by 



(Ml 

M2 

— M^sin6'vi/cos/? M^cos^^/cos/? 

\ M^sin6'vi/sin/3 — M2Cos6'vi/sin/3 



—Mzsm6wcosl3 
M^cos6'vi/ COS/9 




M^sin^VKsin/? \ 
-M^cos6'vi/sin/5 

/ 



(4) 



where we have used the convention of refs.|0, which we will follow throughout. As- 
suming Grand Unification of the gauge couplings implies a relation between the SUSY 
breaking gaugino masses Mi and M2 as well as the gluino mass rriq = |M; 



3|- 



Ml = - tan^6'vi/M2 



5ap 



3as cos^O 



w 



(5) 



where is the electromagnetic coupling constant and is the strong coupling. The 
angle /3 in eq.(^ is defined via tan/? = {H2) / (Hi). Finally, the parameter /i describes the 
supersymmetric contribution to the Higgs(ino) masses. 



In general the LSP is a complicated mixture P3| , ^ P| of the 4 current eigenstates; 
in our numerical calculations we take full account of this mixing by diagonalizing the 
mass matrix numerically. However, in the limit \Mi \ + ^ Mz this diagonalization 
can quite easily be carried out perturbatively. Since this proves helpful for a qualitative 
understanding of our numerical results, we list the eigenvalues rrii and eigenvectors of 
the mass matrix (H), keeping terms up to 0{Mz)- 



mi = Mi; 




MzsinOwicospMi + sin/?/i) Mzsin6'vi/(sin/5Mi + cos/3/i)' 



/i2 - 



Mzcos9w{cospM2 + sin/3/i) M2Cos6'iy(sin/3M2 + cos/5/i)\ 
mg = /i(l + S); 

1 / M^sin^pi/ (cos/3 + sin/3) M^cos^vk (cos/3 + sin/3) \ 

= m;^ ' j^^2 'I'-i + ^j! 

mi = -/i(l + 6'); 

1 f MzsinOwicosP — sinP) M^cos6'vi/(sin/? + cos/3) 



64 



1,1 + e' 



(6a) 



(6b) 



(6c) 



(6d) 



Ml + /i M2 + H 

Note that eq.(|^) implies \m2\ > \mi\] 62 therefore never corresponds to the LSP. In the 
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limit where either \Mi\ or (or both) is much bigger than M^, the LSP is therefore an 
almost pure bino (ci), or an almost pure antisymmetric or symmetric higgsino (e3 4).0 

The corrections 5, 5',e and e' are formally (9(M|); however, they can be numerically 
quite important for a light higgsino-like LSP. They are given by 

^ _ M|(l + sin2/3) ( sm^Ow cos%.\ 
^ M|(l-sin2/?) / sin^^H- cos^\ _ 



. cos2/? 
1 + sin2/?' 



(7c) 



e' = (7d) 
1 - sin2/? ^ ' 

After inclusion of these terms one therefore finds that the mass splitting between the two 
higgsino-like states is given by (in the limit |Mi| ^ |/i|): 

""^^1 - 1"^^" = 3|M,| V + ^^M^' il^V • 
More details about neutralino masses and mixings can be found, e.g. , in refs. |1T7| , ^ |43| , ^ . 



What are the final states into which a pair of neutralinos can annihilate? Here we only 
include 2-body final states that can be produced in leading order of perturbation theory. 
From unsuccessful sparticle searches at LEP as well as the bound rrig > 120 GeV 
that follows from the preliminary gluino search limit of the CDF collaboration 
inclusion of "cascade decays" [0, one can derive |^ the bound > 20 GeV. 
that the annihilation xx ~^ ff is always kinematically allowed for all light SM fermions up 
to and including b quarks. For heavier neutralinos annihilation into a pair of gauge bosons 
also has to be included ||18|, |19[ . In addition the model contains at least one neutral scalar 



4g] after 
We see 



Higgs boson h with mass not much above Mz, even after inclusion of radiative corrections 



30|; the second neutral scalar H, the pseudoscalar P and the charged Higgs boson 



can also be accessible. In general one therefore also has to include annihilation into two 



Higgs bosons |T^, |T^, as well as into one Higgs and one gauge boson. 

We will now discuss annihilation into these final states in a little more detail, assuming 
X to be either a nearly pure bino or a nearly pure higgsino. Here we only give symbolic 
expressions for the matrix elements, which allow to estimate the order of magnitude of 
the various contributions; exact expressions for the cross sections for a general x state 
are listed in Appendix A. Since we expand the annihilation cross section of eq.(|^) only 
up to O(f^), we only have to include s and p wave contributions, s wave contributions 



*Unless \Mi\ ~ in which case strong higgsino-bino mixing occurs. Such mixed states always lead 
to very small relic densities |2|, 

^ We note in passing that this implies a freeze-out temperature Tp > 1 GeV, well above the temperature 
where the quark-hadron phase transition is expected to occur. Our results do therefore not depend on 
the exact value of the critical temperature for this phase transition, in contrast to the case of a very light 

LSP M- 
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start at 0{v^), but also contain terms that contribute to eq.(Q) via interference 

with the 0{v^) terms, p wave matrix elements start at 0{v), so that we only need the 
leading term in the expansion. Of course, there is no interference between s and p wave 
contributions, and hence no 0{v) term in eq.(^. Notice finally that Fermi statistics forces 
the s wave state of two identical Majorana fermions to have CP = —1, while the p wave 
has CP = +1; the same argument also implies that the s wave has to have total angular 
momentum J = 0. 



XX ff 



This reaction proceeds via the s-channel exchange of a Z or Higgs boson, as well as via 
sfermion exchange in the t— channel. Each chirality state of / has its own superpartner 
with in general different mass, which mix |^ if the Yukawa coupling of / is not negli- 
gible. In the following expressions summation over both sfermions is always understood. 
Note that both the Z — f — f and fermion-sfermion-gaugino couplings conserve chirality; 
the sfermion and Z exchange contributions to the s wave matrix element Aig are there- 
fore proportional to the mass of the final state fermions. Contributions from Higgs 
exchange, from the higgsino-sfermion-fermion Yukawa interactions, and from sfermion 
mixing violate chirality, but have an explicit factor of rrif. The coefficient a in the expan- 
sion of the annihilation cross section is therefore always proportional to for this 
final state, independent of the composition of x- Moreover, since the CP quantum num- 
ber of the exchanged Higgs boson must match that of the initial state, only P exchange 
contributes to Aig, while h and H exchange contribute to Aip. Since Aip only contributes 
to the coefficient b in eq.(0), which is suppressed by a factor 3/xp ~ 0.1 —0.2, P exchange 
is potentially much more important than the contribution from the scalar Higgs bosons. 
For a bino-like LSP the matrix elements thus have the structure 



bino 



Cl 



Yf + C2 



Ml 



+C3 



1 



mt 



bino 



OC g''^v 



Ml + fi 4m2 — rrip + impTp 

2 



(9a) 



di 



m: 



+ d2 



Ml 



mr 



Ml - /i2 4m2 - M| + iMzTz 



1=1 



Mi + fiAml-Mfj^+tMnTH^ 



(9b) 



Here the Cj, di are numerical constants of order 1, and g' is the U{1)y gauge coupling. 
C3 and the d^^i describe the Higgs-// couplings, where we have introduced the notation 
|2| ifi = if, H2 = h; in certain cases some of these couplings can be enhanced |0 by a 



factor tan/3 (or suppressed by cot/3). 

We see from eqs.(^) that the s— channel diagrams are all suppressed by small couplings. 
The Z boson couples to neutralinos only via their higgsino components; more precisely. 



9zxiXj ^ *2«,3ej,3 — 64,46^4, 



(10) 
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where Cj^; is the /— th component of Cj, see eqs.(|D. For a bino-hke state this couphng 
is doubly suppressed, as indicated in eqs.(^. The Higgs boson couphngs to neutrahnos 
originate from the Higgs-higgsino-gaugino gauge interaction terms of the unmixed La- 
grangian; for a bino-hke neutrahno this couphng involves therefore only one factor of the 
(small) higgsino component. Finally, the sfermion exchange contributions can only be 
suppressed by choosing m| ^ ml in this case. 



For a higgsino-like LSP, eqs.(|) become: 



niggsino 



Mz „ rrif 



1 



Ml + 4m2 — mp + impTp 



;iia) 



higgsino 



oc ig' + g")v 



+c: 



mr 



Mz, 

2 



m 



X 



fiM Ami -^z + iMzTz 



E4 



X 

mj 



3,i 



i=l 



Ml + fi Ami ~ "^P + impTp 



(lib) 



Notice that eqs. ([TI|) also get contributions from SU{2) gauge interactions, which enter 
eqs.(|^) only in higher orders in Mz/{M + /i). On the other hand, the sfermion exchange 
contribution is now suppressed by either the small gaugino component of the LSP, see 
eqs.(6c),(6d), or by a power of the Yukawa couphng; of course, for / = t the latter is 
hardly a suppression, and can even be an enhancement if the top quark is heavy. In 
contrast, Higgs exchange contributes at the same order (in j^^^^ ) to the annihilation of 
bino-like and higgsino-like LSPs.i Finally, the Z exchange contribution again behaves 
quite differently in the two cases of eqs.(P) and (|ll|): While in the former case, it decreases 
quadratically with the mass of the heavy neutralinos (whose mass is oc /i, in this case), 
eqs.(0) and ( p!Ol) show that the Zxx coupling for a higgsino-like LSP decreases only 
linearly with the mass of the heavy neutralinos (with mass oc M)m Moreover, eqs.(|^) and 
(|10|) imply that the off-diagonal Zxx' coupling is not suppressed at all if x = es, x' = ^4 
or vice versa. If Z exchange gives the dominant contribution (which is true for light 
higgsinos, except in the vicinity of the h pole), one therefore has 



(^ixx' ff) oc 



(^{xx ff); 



(12) 



^The fact that the Higgs couphng to a mixed neutrahno is unsuppressed explains to a large part the 
big annihilation cross sections, and hence small relic densities, of this kind of LSP. 

§Unitarity implies that the Zxx coupling must be suppressed at least like as — > oo. The Z 

exchange contribution to Aig behaves like g^^^ra/TO^/M^. This contribution cannot be cancelled by / 
or P exchange because one can always choose m^,mp ^ to^; since this does not increase any couplings, 
the heavy states simply decouple in this limit. The Z coupling itself therefore has to compensate for the 
factor of m^. 
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In this case xx' co-annihilation [37| can be quite important |26|; we will come back to 
this point later. 



XX ^ w^w-, zz 

These final states can be produced via t— channel chargino or neutralino exchange, as well 
as s— channel exchange of scalar Higgs bosons; in case of the W^W~ final state s— channel 
Z exchange also contributes. It is important to realize here that the cross sections behave 
quite differently for longitudinal and transverse gauge bosons. For each longitudinal gauge 
boson the amplitude gets an enhancement factor 7y = EyjMy ~ m^jMy. Unitarity then 
requires strong cancellations between different contributions to the matrix element if the 
couplings of X to are not suppressed. On the other hand, these enhancement factors 
can give finite matrix elements in the limit oo even if the couplings do vanish in this 



limit. These effects can also be understood from the equivalence theorem [|50[, which states 
that in the high energy limit, the matrix element for the production of a longitudinal gauge 
boson is identical to the one for the production of the would-be Goldstone boson that gets 
"eaten" when the gauge boson acquires its mass. The couplings of these Goldstone modes 
to neutralinos and charginos originate from the Higgs-higgsino-gaugino gauge interactions 
of the unmixed Lagrangian. This means that a pair of would-be Goldstone bosons can 
be produced with full gauge strength from pure gaugino as well as pure higgsino initial 
states, by exchange of higgsinos or gauginos, respectively. On the other hand, since the 
relevant couplings are gauge couplings, the matrix element must be well behaved in the 
limit where any mass becomes very large. Notice finally that neither a pair of longitudinal 
gauge bosons nor a combination of one longitudinal and one transverse gauge bsoson can 
exist in a J = state with CP = —1; these final states are therefore only accessible to 
the p— wave initial state. 

For a bino-like LSP, the coupling to gauge bosons does indeed decrease like l/m^, see 
eq.(6a). The above discussion then shows that only the amplitude for the production of 
two longitudinal gauge bosons survives, which is purely wave: {V = W, Z): 



M| Mz m^My 



2 



(13) 

Because of the enhancement factor 7^, the annihilation of heavy bino-like LSPs into 
(longitudinal) gauge bosons does not vanish for 00, unless one has \Mi\ -C 

Considering |]19[ an exact bino state therefore does not give the right answer in this 
case. As discussed above, this can also be understood from the equivalence theorem, 
since the production of two Goldstone modes can only be suppressed by making the 
exchanged higgsino very heavy. (This argument is made more rigorous in Appendix 
B, where numerical factors from the diagonalization of the Higgs sector etc. are treated 
properly.) Finally, we mention that the coefficient dr, i is often quite small, since the heavy 



Higgs scalar decouples |^ from W and Z bosons when its mass is large. In contrast the 
exchange of the light Higgs boson actually gives the dominant contribution in the limit 
> |Mi|. 
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For a higgsino-like LSP one has 



MsixX -V^V^)kiggsino « (r + 9"K, (14a) 



Mz m^Mv ml 

^'^M + /i 4m2 - ml,^ + ^m^TH, 



(14b) 



Notice that in this case there is no propagator suppression of the t— channel diagrams for 
VtVt production, because the exchanged (other) higgsino state is almost mass degenerate 
with the LSP. Since the (off-diagonal) couplings to gauge bosons are not suppressed in 
this case, cancellations between t— channel (and Z exchange, for V = W) diagrams are 
necessary to restore unitarity for V^Vl production. The equivalence theorem shows that in 
the limit \Mi\ ^ the contribution from longitudinal gauge bosons is even suppressed, 
since the production of Goldstone bosons necessitates the exchange of heavy gaugino-like 
states, or the exchange of Z or Higgs bosons whose diagonal couplings are also suppressed 
in this limit. Moreover, the production of VlVt final states is suppressed by neutralino 
mixing factors. However, the production of transverse gauge bosons, and hence the total 
cross section, is not suppressed in the limit where the masses of the heavier neutralinos 
becomes very large; this is again in contrast to the case of a bino-like LSP.0 

XX Zh 

This final state can be produced via neutralino exchange in the t— channel as well as 
s— channel exchange of Z and P bosons. Notice that now a longitudinal Z boson can be 
produced both from the s— and p— wave; since we always need at least one neutralino 
mixing factor, the matrix element would go to zero in the limit of large LSP mass without 
the enhancement factor 7^:. For the case of a bino-like LSP, the matrix elements have 
the form: 



Ms{xX^Zh)\,,^^^g- 



1 

Cq + C7- ^ 



4m2 — Trip + iiripTp 



(15a) 



9 

M.ixx -^/^)Li„o ^ 9"vd,j^^^^. (15b) 

Both cg and get contributions from neutralino as well as Z exchange diagrams; in the 
important limit mp S> M|, is strongly suppressed, due to a cancellation between 
the two classes of diagrams. In this limit the ZhP coupling also becomes small 
suppressing the P exchange contribution as indicated. Only the p— wave amplitude (15b) 
therefore survives in the limit mp ^ M|. Notice, however, that the amplitude as a whole 
does not vanish in the limit of large sparticle masses and small neutralino mixing, unless 



^The VlVl final states have not been treated properly in ref.[|9|. However, since they only contribute 
to p— wave annihilation, they are numerically not very important. This is even true for bino-like LSP's 
where the production of transverse gauge bosons does not contribute, since here the total annihilation 
cross section is dominated by the // final state. 
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IfJ'l ^ \Mi\; once again this is due to the production of a longitudinal gauge boson, giving 
rise to an enhancement factor m^/Mz- 

For a higgsino-like LSP eqs. ([l5|) become: 



HI 



M + fi 



rr?p + M| 4m?, — mp + iTpnip 



(I6a) 



2 

m: 



M.ixX -^/^)Liggsino « (^?' + 9>d7jpf^- (16b) 

Note that in this case the 0{v^) term from the t— channel and Z exchange diagrams is 
not suppressed for mp 3> M|. Just as in the case of a bino-like LSP the total amplitude 
is only suppressed if the heavier neutralinos are much heavier than m^, i.e. if \Mi\ ^ 
in this case. 

XX hh 

Here only t— channel neutralino exchange and s— channel scalar Higgs exchange diagrams 
contribute. Since two identical scalars cannot be in a state with J = and CP = — 1, 
annihilation can only proceed from the p— wave. The amplitude thus has the general form: 



Mpixx -^hh) (xg V 



, , Ml Mz Mzm^ 



'M + fi "M2-/i2 ^ '"''M + fiAml-ml +tmHTH, 

(17) 

This form holds for both bino-like and higgsino-like LSP's. The contribution oc dg comes 
from the exchange of the heavier neutralinos, which occurs with full gauge strength but 
is suppressed by small propagators; the term oc dg originates from neutralino mixing. In 
case of a bino-like LSP the coefficient ds is suppressed if tan/? ^ 1 and ^ \Mi\: 

Ml + /X sin2/? 



This possible additional suppression is absent for the case of a higgsino-like LSP; in this 
case the amplitude also gets contributions from SU{2) gauge interactions, as do all other 
higgsino annihilation amplitudes. 



This concludes our qualitative discussion of the annihilation matrix elements for the most 
important final states. In principle the heavier Higgs bosons H, P, could also be 
produced [^, |19| in xx annihilation, either in pairs or in association with a gauge boson. 
However, we will see in the next subsection that in minimal supergravity models these 
heavy states are usually not accessible. We do therefore not discuss the relevant matrix 
elements here; of course, they are included in the list of matrix elements in Appendix A. 
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2c. The particle spectrum 

The basic assumption of minimal Supergravity (SUGRA) models is that supersymmetry 
breaking can be described 0] in terms of just three parameters: A universal scalar mass 
m, a universal gaugino mass M, and a universal parameter A characterizing the strength 
of nonsupersymmetric trilinear scalar interactions. If the particle spectrum is restricted to 
that of the MSSM, which we always assume in this paper, one in addition has to introduce 
a supersymmetric contribution /i to Higgs boson and higgsino masses; the masses /ii and 
/i2 of the two Higgs doublets are then given by 

2 2 2,2 /inN 

jj^ = IJ,2 = m + n , (19) 

while the Higgs mixing term /x| is given by 

/i^ = (A - m) ■ m. (20) 

This simple form of the particle spectrum is assumed to emerge after integrating out the 
fields of the "hidden sector" where local supersymmetry is broken spontaneously. Since 
the decoupling of these fields occurs at the Planck or GUT scale the spectrum will be this 
simple only at ultra-high energies Q > Mx', in particular, eq.([T9|) will only hold at these 
very high energies. 

Of course, present day experiments, as well as DM annihilation, occur at much smaller 
energy scales. The particle spectrum at energies of the order of the weak or sparticle 
mass scale can be obtained by solving a set of coupled renormalization group equations 
(RGE)[0; for given m, M, A, /i and Yukawa couplings ht, hb (with hb = hr at the GUT 
scale), minimal SUGRA specifies the boundary conditions at Q = Mx, which uniquely 
determine the particle spectrum at lower energies. In this scheme leading logarithms 
are automatically summed, i.e. all terms of order (a/vr In Mx/M^)" are automatically 
included, where a is a generic gauge or Yukawa coupling. It has been recognized quite 
early |T0| , [5T[1 that the radiative corrections described by the RGE can induce spontaneous 



breaking of the SU{2) x f/(l)y gauge symmetry by driving some combination of the 
squared Higgs mass parameters /i^ of eq. ( pJ]) to negative values. This is due to the effect 
of the Yukawa couplings, which tend to reduce the squared masses of scalar fields. The 
Yukawa sector therefore plays a crucial role in these models. 

In a recent paper we studied radiative gauge symmetry breaking and the particle 
spectrum in minimal SUGRA in some detail, taking care to incorporate the effects of 
the Yukawa couplings of the h quark and r lepton, which can be quite important if 
|tan/3| ^ 1. We later showed |Q how to incorporate the "finite" (i.e., without InMx/M^ 
enhancement) radiative corrections to the Higgs sector in this scheme. In particular we 
demonstrated that for most purposes these radiative corrections can be made negligibly 
small by a proper choice of the scale Qo where the RG running is terminated; the only 
exception is the mass mh of the light Higgs scalar where the corrections have to be included 
explicitly. (A similar result had been obtained previously in ref. |53|].) For more details we 



refer the reader to these papers, as well as to earlier work on this subject ||Ty, ^; here 
we only give a brief summary of the relevant properties of the spectrum. 

As stated above, the squared mass of all sfermions at the GUT scale is simply given 
by rn? . In case of the superpartners of the first two generations the only sizable radiative 
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corrections to the masses involve the gauge interactions. The breakdown on SU{2) x f/(l)y 
gauge symmetry also has some effect on sfermion masses. All these contributions can quite 
easily be computed analytically; one has 



mj^ =m^ + diM^ + cos2/3(4, - Qi sin^^vK)M|, (21) 

where /s^j and Qi are the third component of the weak isospin and electric charge of the 
sfermion /j, respectively. The positive constant di is determined by the gauge quantum 
numbers of the sfermion; numerically, it is about 6 for squarks, 0.5 for SU{2) doublet 
sleptons and 0.15 for SU{2) singlet sleptons with hypercharge 1. 

The gaugino masses are all equal to M at scale Mx; their Q dependence is identical 
to that of the gauge couplings a^: 

M.(Q) = ^M, (22) 

which immediately implies eq.(|^). Numerically, M3 ~ 3M, M2 — 0.84M and Mi ~ 
0.43M.E 

Eqs.(pl|) and (^) have been taken into account in some previous analyses |22 



^ of LSP relic densities. Unfortunately the effects of the Yukawa couplings are 
not so easily treated analytically. These effects are important for the masses of third 
generation sfermions as well as Higgs bosons. Indeed, Yukawa couplings affect sfermion 
masses already at tree level |^; they lead to mixing between SU{2) doublet and singlet 



sfermions, and give rise to additional supersymmetric diagonal mass terms. These effects 
are especially important for stop squarks, which obviously have the largest Yukawa cou- 
plings; this has been included in the analysis of ref. EH]. However, for |tan/5| ^ 1, sbottom 



and especially stau mixing also becomes important . Moreover, the Yukawa couplings 
reduce the nonsupersymmetric diagonal mass terms from the values predicted by eq.(pl|). 
The net effect is that the lighter stop and stau eigenstates can be substantially lighter 
than the other squarks and sleptons, respectively; indeed, no strict lower bound on these 
masses could be given even if the masses of first generation sfermions were known. The 
maximal reduction of the light sbottom mass is not quite as large, but can still amount 
to 20 - 30 % if |tan/3| ~ mt/rrib- Since tree-level and loop effects tend to cancel for 
the heavier eigenstates of third generation sfermions, their masses are usually not very 
different from those of the corresponding sfermions of the first two generations. 

The Yukawa couplings also affect the Higgs masses via the RGE; as explained above, 
this effect is at the heart of the radiative gauge symmetry breaking mechanism. First 
of all, it should be noted that the masses of the heavier Higgs bosons are related to the 
sfermion masses via eq. ([l9|) ; this equation only holds at scale Mx, but it shows that in 
SUGRA models one cannot treat Higgs boson and sfermion masses as independent free 
parameters. As well known ||5^, the masses of the P, H and states are essentially 
determined by mp, which is simply given by 

ml = filiQo) + fiUQo); (23) 



*In our numerical calculations we use proper on-shell masses, i.e. mj = ni^ [Q — rrij ) and similar for 
gluinos. The coefBcients in eqs.(^l|) and (^) then depend on the masses themselves, rather than being 
simple constants. This can change the masses of strongly interacting sparticles by as much as 20%; e.g. 
for 1 TeV gluinos, rrig ~ 2.5|M|, rather than 3|M|. 
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as discussed above, eq. (p3D also holds to good approximation after inclusion of 1-loop 
radiative corrections if the scale Qq where the RG running is terminated is chosen properly, 
i.e. Qo ~ rrig. Notice that only the Higgs doublet H2 couples to top quarks; moreover, 
the requirement 

{H",f + {HlY = 2M^/g' (24) 



immediately implies flOl fJ^KQo) > —Mz/2. If the b and r Yukawa couplings are neglected, 
the nonsupersymmetric contribution to /if runs just like an SU{2) doublet slepton mass, 
i.e. it increases when going from Q = Mx to Q = Qq; in addition the positive super- 
symmetric contribution /.^ has to be added 0. This argument shows that mp can only 
be smaller than the slepton masses if the b and r Yukawa couplings are sizable, i.e. if 
I tan/3 1 ^ 1. The exact numerical expression for mp can be approximated by K 



ml = ^ (cot/5 - 1) + (m^ + 0.52M2 + f^\Qo)) ^ - — ^ j , (25) 

where a ~ 1/45 — 1/35 depends on the ratios M/m^ A/m as well as on the top mass mj.i 
The mass of the lightest Higgs scalar h in general depends on all parameters of the 
model in a complicated way. However, in the important limit mp ^ M| the situation 
simplifies greatly, and one finds [p^ : 



ml = Ml cos' 2/5 + A22 sin'p + o(^,^]. (26) 

y mf, mf J 

Here A22 describes the leading radiative corrections from top-stop loops. It grows like the 
fourth power of mt, but depends only rather mildly on the values of the SUSY breaking 
parameters; in the limit m? ^ mf it grows oc In mq/mt. 

The SUGRA-imposed constraint that is most difficult to treat analytically follows from 
the almost obvious observation that in the radiative gauge symmetry breaking scenario 
the vevs of the Higgs fields can be computed from the input parameters at the GUT 
scale, i.e. (-^12) functions of m, M, A, /x and the set of Yukawa couplings ht, hb, h^-. 
The condition ( p4D therefore leads to a relation between these parameters. This relation 
cannot be expressed in closed form once the effects from r are included, but it can 
approximately be written as 



'"'^^"^ - tan'^-1 ^^ ~ " " ^^'^^ 

where X2 describes the effect of the RG-running due to the top quark Yukawa coupling. 
An approximate expression for X2 isi 

\^ ] n n™2 , o '7A/r2 



150 GeV. 



) I 0.9m2 + 2.7M2 + 1 - (^—^^^ (0.24^^ + MA) (28 



^^Note that this is the square of the real parameter /i, which renormahzes multiphcatively; /i^ is 
therefore indeed always positive. 

^In ref. we gave a somewhat larger numerical value of a, because we under-estimated the effect 
of the running of the b quark mass between Qq ^^nd mf,. Since a smaller to;,((3o) sXso implies a larger 
upper bound on |tan/?|, our results for sparticle masses, including mj,, remain valid if one rescales tan/3 
1.15 tan/? in the region | tan/3 1 ^ 1. 

^In ref.[^ we gave a somewhat larger expression for X2, since we had chosen a rather small value for 
go, i.e. Qo = Mz] eq.(|2|) is valid for Qo ^ 300 GeV. 
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Eq.(l^) usually works to 10% accuracy, but eq.(^) might deviate by as much as 20% 
from the exact numerical result; nevertheless these expressions are quite useful to gain 
some insight into the relation between the input parameters. 

Strictly speaking, eq.(|27|) does not provide the searched-for relation between input 
parameters at scale Mx on the one hand and Mz on the other, since it still depends 
on tan/3, which is itself a complicated function of the input parameters. (In fact, this 
is where the dominant effect from the other Yukawa couplings enters.) However, /i^ 
becomes essentially independent of tan/3 if tan^/3 3> 1, in practice for |tan/3| > 3 or so. 
Another complication arises because a given set of m, M, A and often allows up to 



three different solutions ||3^ of the equation that determines tan/3 (and hence the Yukawa 
couplings); these solutions differ in both sign and magnitude of tan/3.0 We therefore 
present results always for fixed tan/3, M and m, rather than fixing the values of the 
SUSY breaking parameters. Since consistent solutions only exist for ht > hf,, the allowed 
range of tan/3 in this model is restricted to 1 < |tan/3| < mt/mb. 

Unfortunately another two-fold ambiguity occurs when /i and A are adjusted such 
that Mz and tan/? have their desired values. In one of these solutions /i is of the order 
of m and M. However, unless is very large there is also a second solution with 
/i^ <^ m^, M^; indeed, this is the "small /i" solution which has been discussed in the 



first analyses of radiative gauge symmetry breaking. However, by now this kind of 



solution is quite severely constrained. First of all, we know from SUSY searches at LEP 



45| that |/i((5o)| > 40 GeV; since these solutions typically have much smaller than 
m and M, this constraint implies that most sparticles must be quite heavy for these 
solutions to be acceptable. Moreover, for nit > 155 GeV the small—// solutions disappear 
altogether, since then the effect of the top Yukawa coupling always drives /i^Qo) below 
— M|/2 unless it receives a sizable, positive contribution see eq.(|9]). 

At this point some comments on fine tuning might be appropriate. It should be quite 
obvious that the "natural" scale for the vevs of the Higgs fields is set by the dimensionful 
parameters of the Higgs potential, which in turn are roughly of the order of typical 
sparticle masses, as can be seen from eqs. ([23|) and (0). Therefore some fine tuning 
of parameters will be necessary |16| to achieve M| -C + M^. This provides a 



strong argument that sparticles should not be much heavier than 0{1 TeV), but this 
argument cannot be translated into strict upper bounds on sparticle masses. Since one 
of the motivations of this study is to see whether cosmology might provide us with such 
bounds, we do not impose any "a priori" upper bounds on the SUSY breaking parameters 
in our analysis. One might also argue that large ratios of the mass parameters of the 
model are "unnatural", but the Yukawa sector shows that large ratios of "fundamental" 
parameters can indeed occur. Of course, one ultimately hopes to understand the origin of 
the dimensionful parameters m, M, A and fi better, e.g. in the framework of superstring 



theories |£5[. However, at present it seems safer to pursue an "agnostic" approach, and 



try to cover the entire experimentally allowed parameter space. 



^In this scheme the sign of tan/3 is is determined dynamically via the RGE. On the other hand, the 
couplings listed in ref.Q] and Q have been derived under the assumption tan/3 > 0. Fortunately the 
masses and mixings of the neutralino and chargino eigenstates only depend on the sign of the product 
M ■ n ■ tan/3; a sign in tan/3 can therefore always be "rotated" into M or fi. 
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This allowed region is defined via the following experimental constraints: 



m- - - > 45 GeV; (29a) 

> 40 GeV; (29b) 

^ Br{Z ^ x-X°) < 5 ■ 10-^; (29c) 

T{Z ^ XiXi) < 12 MeV; (29d) 

m-g > 120 GeV; (29e) 

rrif-^ > m^. (29f) 



The bounds (29a)-(29d) directly follow from LEP limits on sparticle production |^ 



as 



well as on the invisible width of the Z boson; here and stand for a generic chargino 
and neutralino state, with x5 = X- The bound (29e) is a conservative interpretation of the 
preliminary CDF search limits |^ after inclusion of cascade decays [0. Finally, (29f) 



follows directly from the requirement that the LSP should not be charged, as discussed in 
the introduction. Further experimental constraints follow from the unsuccessful search of 
Higgs bosons; we have incorporated a parametrization of the ALEPH bound |^ in our 
list of conditions. 

In addition to imposing these experimental constraints, we also discard combinations 
of parameters that lead to deeper lying minima of the scalar potential that break charge 
and/or color; this requirement excludes combinations with / {m? + M"^) ^ 1 p7 



Finally, we demand that the scalar potential should be bounded from below at scale 
Q = Mx, which implies + fil > 2|/i3|; eq.(pOD shows that this excludes the region 




l-Bl / B'^ \B\ 

y - - K HMx)\ < y + ^/^ - 1, (30) 

if \B\ = |yl — m| > 2m. This constraint is effective at small |tan/3|, which implies large /i, 
see eq. (pT]) , and thus large and positive squared Higgs mass parameters (|2T|). One then 
needs large |A/m| to achieve spontaneous gauge symmetry breaking, since this accelerates 
the RG running of the mass parameters; however, large \B/m\ also imply a large excluded 
region (|30|). 



3. Examples 

We are now in a position to present some numerical results. The discussion of sec. 2c 
showed that the model has 4 free parameters, which we chose to be m, M, rrit and tan/3; /j, 
and A are then fixed by the equations describing the minimization of the Higgs potential, 
up to a possible discrete ambiguity. Without loss of generality rrit and m can be chosen 
to be positive, but M and tan/5 can have either sign. 

Figs. 1 a-d show a first partial exploration of the parameter space of the model. In 
these figures we have fixed m = 300 GeV, which leads to cosmologically interesting DM 
densities for a wide range of the remaining parameters. In addition, in each figure we have 
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kept rrit and tan/3 fixed, and varied M. We find that for the chosen values of parameters 
only one experimentally allowed solution for A and fi exists. Moreover, since in all these 
cases the mass of the LSP increases monotonically with |M|, we present our results as a 
function of m^; this simplifies the identification of the various s-channel poles and of the 
thresholds where new annihilation channels open up. The starting point of all curves in 
figs. 1 is determined by the LEP constraints (29a), (29c) and (29d). 

We see that in the limit of large all curves become almost identical. The reason 
is that in this case we always have > \Mi\ ^ Mz, so that the LSP is bino-like. The 
dominant annihilation channel is then XX ~^ via I exchange, where / is any lepton. 
Eq.(l) shows that squark exchange is suppressed 



2| by their large mass (^), and the Z 
and Higgs exchange contributions are suppressed by small couplings. As discussed in sec. 
2b, the production of longitudinal gauge bosons Vl and of the light Higgs boson h are not 
suppressed by powers of Mz/m^, but are suppressed by powers of |Mi//i| ~ 1/4 — 1/2 
for the examples of figs. 1; moreover, all neutralino couplings relevant for these final 
states originate (either directly or via the equivalence theorem) from the Higgs-higgsino- 
gaugino U{1)y gauge interaction, which involve fields with hypercharge \Y\ = 1/2, while 
the SU{2) singlet leptons have Y = 1. The gauge and Higgs boson final states therefore 
only contribute a few % in the region > 200 GeV (> 300 GeV for the dashed curve 
in fig. Ic; see below). 

The behaviour of the curves in this region can therefore be understood semi-quantitatively 
from the Ir exchange contribution alone: 



oc 



m: 



rrij + 



+ 



{rrij + m?,)^ 



(31) 



where the parenthesis results from the Taylor expansion of the propagator in powers of 
the velocity v. SUGRA predicts m3 ~ + 0.83m^ for the given case of a bino-like LSP, 
see eqs.(pT|) and (p2D; eq. ([3T|) then leads to a maximum of the annihilation cross section, 
i.e. a minimum of the relic density, at ~ 0.6m ~ 180 GeV for the parameters of figs. 
1. This maximum is very broad; the Ipt exchange contribution falls to 50% of its maximal 
value at ^ 1.52m, just beyond the end of the region shown in figs. 1. Indeed, in 
Figs. la,b, flh'^ at m^ = 430 GeV is about twice as large as at the minimum. Eq.(PTD 
also predicts the annihilation cross section to fall at small values of m^, dropping to half 
the maximum value at ~ 0.25m. However, in many cases our assumption that Ir 
exchange dominates the total annihilation cross section is no longer valid in this region. 

Going towards smaller values of m^, the first prominent structure one encounters is the 
tt threshold. It is most prominent for the dashed curve in fig. Ic, since this combination 
of parameters leads to the smallest value of for a given value of M, a smaller 
means larger higgsino admixtures to the LSP, see eq.(6a), and hence larger couplings to 
h and Z bosons. For most of cases shown in figs. 1 the contribution of the ti final state 
is suppressed by destructive interference between t and Z exchange contributions, which 
have approximately equal magnitude but opposite signs here. On the other hand, top 
production is not p— wave suppressed unless 3> mf, unlike the production of massless 
fermions. The importance of the tt final state compared to light fermions is therefore 
enhanced by a relative factor xf/3 ~ 10, see eq.(|^). 
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In the region below the ti threshold the differences between the various curves start 
to become more pronounced. Many of these differences can be understood from eqs. (p7D 
and (0), which show that increasing rrit and decreasing tan/3 both imply larger values of 
/i, which leads to smaller couplings of the LSP to Higgs and gauge bosons, and suppresses 
contributions from the t— channel exchange of higgsino-like, heavier neutralinos. This 
explains why the WW and ZZ thresholds, as well as the minima at = Mz/2 and at 

= m/i/2, are more prominent in fig. la than in lb, and why Qh"^ at small m^, where 
Z exchange diagrams dominate the annihilation cross section, is considerably smaller in 
figs. lc,d than in la,b. Moreover, rrih increases with increasing rrit and increasing |tan/3|, 
so that the position of the /i— pole tends to move to larger values of as we go from fig. 
la to Id. 

The depth and width of the minimum at = mh/2 depends quite sensitively on the 
choice of parameters. Eq.(P3|) shows that in figs. 1 we always have mp 3> M|; in this 
limit the hxx coupling |^ becomes for a bino-like LSP: 



_ g' Mzsm9w{M,+fism2P) ( M| M|\ 

Notice that the two terms in the numerator tend to cancel if Mi ■ fi ■ tan/3 < 0; this 
explains why the minimum at = mh/2 is narrower and shallower for the dashed 
curves in figs. 1 than for the solid ones. However, we remind the reader that estimating 
Qh"^ from eqs.(|T])-(^ can lead to large errors in the vicinity of a very narrow pole. A 



more careful treatment would lead to shallower minima, which are broadened in the 
region below mh/2] if mh > Mz, the relative maximum between the two minima should 
therefore also be somewhat lower than indicated in figs. 1. On the other hand, the h— 
pole clearly affects only a very limited region of parameter space; our overall conclusions 
do therefore not depend on an accurate treatment of this pole. 

Finally, the strength of the Zxx coupling also depends quite sensitively on the choice 
of parameters, including their signs. In this case the ordering of the curves for M > 



and M < even depends on tan/3. For tan/? = 2, a cancellation occurs Q in the 
neutralino (and chargino) mass matrix if Mi ■ ^ ■ tan/3 > 0; the LSP then has substantial 
higgsino components if m^ < Mz/2. In contrast, for small tan/3 and Mi ■ fi ■ tan/3 < 
the LSP remains dominantly a gaugino even if is very small. On the other hand, 
for M < 0, tan/3 = 15, fi{Qo) is quite small due to a strong cancellation in the r.h.s. 



of eq. (|27|) , especially for = 130 GeV, fig. Ic. At larger |M|, i.e. larger m^, this 
cancellation is less complete, but for this choice of parameters x remains dominantly a 
higgsino for m^ < 60 GeV, and reaches 90% bino-content only for > 120 GeV; for 
M > and the same values of mt and tan/3, the LSP has already 97% bino-content at 
this mass. This explains why the two curves in fig. Ic differ quite strongly even at rather 
large values of m^. 

The strong dependence of some of the annihilation cross sections on the model pa- 
rameters is further illustrated by fig. 2, which shows a blow-up (on a linear scale) of 
the VV and Zh threshold region for two of the curves of figs. 1. Both curves show 
small shoulders at the WW and ZZ thresholds. As explained above, these thresholds 
are somewhat less pronounced for the case mt = 160 GeV, due to the larger value of fi. 
However, since for fixed m^, /i only increases by approximately 25% as mt is increased 
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from 130 to 160 GeV, this effect is not very large; see also eq.(|^). This rather small 
change of parameters suffices, however, to reduce the XX ~^ hh cross section by as much 
as a factor of 6! For rrit = 130 GeV, the exchange of the lighter, gaugino-like and heavier, 
higgsino-like neutralinos gives contributions of approximately equal size and equal sign 
to the matrix element. In other words, the terms oc ds and oc dg in eg . {\17\} have about 
equal magnitude here; notice that the larger SU{2) gauge coupling can only enter via dg 
in case of a bino-like LSP. Going to larger rrit does not only increase fi for fixed M, it also 
increases rrih via the radiative correction A22 of eq.(|26|); this necessitates an increase of 
M, and thus a further increase of fi, in order to reach the hh threshold. The contributions 
oc dg are therefore almost negligible for the dashed curve in fig.2. Indeed, the decrease of 
this curve after the maximum at = 62 GeV has nothing to do with any thresholds; 
rather it is caused by the increase of the slepton exchange contribution to the final 
state, see eq.(^). Finally, we remind the reader that a cancellation occurs in the contri- 
bution to the hh final state from higgsino exchange if Mi ■ fi ■ tan/5 < 0, see eq.(|l^); this 
contribution is also suppressed for large tan/? if ^ \Mi\. The hh threshold is therefore 
all but invisible in most curves in figs. 1. 

It has been noted in ref. |^ that our estimate of Qh^ eqs.(|l|)-(|^), becomes unreliable in 
the vicinity of a threshold for a final state which quickly dominates the total annihilation 
cross section. The example given there was exactly the hh final state, which for certain 
combinations of parameters can dominate the total cross section by a large factor. How- 
ever, we find that such a situation never occurs for bino-like or mixed LSP's in minimal 
SUGRA; in this case the hh threshold, as well as all other thresholds, is never much more 
pronounced than for the solid curve in fig. 2, and usually the thresholds are much less 
important, as can be seen from figs. 1. We do therefore not expect the error introduced 
by our approximate treatment to exceed 10% just below threshold, and it should be much 
smaller everywhere else. 

So far we have concentrated on examples where the LSP is dominantly a bino, the 
exception being the dashed curve in fig. Ic. Since the LSP will only be higgsino-like if 

< \Mi\ ~ 0.43|M|, one obviously needs quite large values of \M\ to get a higgsino-like 
LSP with mass substantially above Mz- The SUGRA constraints then imply that one 
also needs large m, since for ^ eqs.(^) and p8| ) always yield > \Mi\ for 
experimentally allowed values of rut. In fig. 3 we have therefore chosen m = 2 TeV, and 
present results for two different choices of M and rrit. In this figure, /j,, A and tan/? all vary 
along the x— axis, with ~ As in figs. 1, we see that Qh^ is essentially independent 
of most parameters if the LSP is a heavy, almost pure state, here an almost pure higgsino. 
Of course, once |yu| > \Mi\ the LSP will become bino-like again, and will thus have a very 
small annihilation cross section, since the very large m implies very heavy sfermions; this 
explains the steep rise of the dashed curve at ~ 430 GeV. 

Notice that this curve does not have a relative minimum in the region ~ \Mi\, 
even though here both the higgsino and gaugino components of the LSP are large; this 
seems to be in conflict with results of refs. [1^, However, in our case eq.(^) implies 
that all Higgs bosons except h are very heavy, so that most final states containing Higgs 
bosons are not accessible; moreover, contributions from P— exchange in the s— channel 
are suppressed, even though the Pxx couphng is large for a mixed LSP. The xX ~^ f f 
cross section does show a maximum in the region where x is a mixed state, due to the 
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contributions from Z and h exchange, but this is not sufficient to compensate the rapid 
decrease of the W and Zh cross sections. 

For < M^y the annihilation of higgsino-hke LSP's is dominated by Z— exchange 
diagrams. Eqs.(|TO|) and (|^ show that gzxx °^ ^^is case; the xx annihilation 

cross section is therefore about 4 times smaller for M = 2 TeV than for M = —1 TeV. 
However, the results of fig. 3 are quite misleading in this region. We had already men- 
tioned above that eqs.(||)-(|^) are not valid close to a threshold where a new channel 



opens up which quickly dominates the total annihilation cross section, as is the case for 
the WW threshold here; this is because we ignored the possibility that LSP's with mass 
below Mw can annihilate into WW pairs if they have sufficient kinetic energy. Therefore 
we have overestimated the relic density in the region just below and at the WW threshold. 

Moreover, as already discussed in sec. 2b, a light higgsino-like LSP always implies 
the existence of a second higgsino-like state x' whose mass is quite close to that of the 
LSP. Since the Zxx' coupling is not suppressed, unlike the Zxx coupling, the xx' co- 
annihilation cross section is much larger than crannixx)^ see eq. ([T^) . Furthermore, xx' 
co-annihilation is not p—wscve suppressed even if the final state fermions are massless. 
Using the formalism of ref. [^, we estimate that inclusion of xx' co-annihilation would 
reduce flh^ hj a. factor 

^ (l + yI^-"^) / (l + e-^'^T- (33) 

Here A = KImsl — |m4|)//x| is given by eq.(H) and e by eq.(^. The exponential factor 
describes the Boltzmann suppression of the x' density at freeze-out; we find Xi;' ~ 25 in 
this case. The factor of Xp/S in front of the second term in the numerator of eq.(P3D 
has been included to estimate the s— wave enhancement (or, more accurately, lack of 
jo— wave suppression) of the co-annihilation process, and the 2 is a statistics factor P7| . 
Numerically we find /sT ~ 4 (80) for |M| = 1 TeV and = 50 (80) GeV; for |M| = 2 TeV 
the corresponding numbers are 230 and 1500, respectively. Of course, these estimates 
could easily be off by a factor of 2 or so; nevertheless, taken together with sub-threshold 
XX annihilation into W pairs, these large suppression factors allow us to conclude that the 
relic density of higgsino-like LSP's will always be uninterestingly small unless > 500 
GeV or so.0 

We had seen in figs. 1 that over a wide region of parameter space the relic density of a 
heavy, bino-like LSP depends only very little on m^, tan/? and the sign of M. However, as 
already pointed out in ref. [^, this is no longer true for very large values of |tan/3|. This is 
illustrated in fig. 4, where we show Qh'^ as a function of tan/3 for fixed m, M and mt] the 
parameters are chosen such that x is bino-like. The full line shows the SUGRA prediction 
including the contributions from the b and r Yukawa couplings to the neutralino-fermion- 
sfermion interactions as well as to the RGE. These latter contributions reduce mp, as 



*0f course, co-annihilation will also occur for > Mw; indeed, the relative mass splitting A 
between the higgsino states, and hence the Boltzmann suppression of the co-annihilation contribution, 
will be (much) smaller than for light higgsinos. On the other hand, (Tann(xx') should not be much bigger 
than crann(xx) ill this case, since crann(xx) is dominated by annihilation into a pair of gauge bosons, which 
occurs with full gauge strength; in the region m^^ > M\y co-annihilation should therefore not change the 
result of fig. 3 much. 
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shown in eg . (pSD ; they also reduce the masses of the hghter b and f eigenstates by mixing 
between SU{2) singlet and doublet states, and by reducing the diagonal entries of their 
mass matrices. 

Fig. 4 shows that both these effects are quite important. When tan/5 is increased 
sufficiently, one eventually has irtp = 2m^; for the parameters of fig. 4 this happens 
at tan/3 ~ 35. This results in a very strong enhancement of the XX ~^ bb, t~^t~ cross 
sections via the exchange of a pseudoscalar Higgs boson. Note that the Pbb and Pt^t~ 
couplings increase cx tan/?, so that the total decay width Fp oc tan^/?. For tan/5 = 35 the 
mass to width ratio of P is therefore similar to that of the Z boson, so that our estimate 
of flh"^ should be quite reliable even in the pole region ||37| . 

The long dashed curve has been obtained by artificially keeping mp constant at the 
value SUGRA predicts for tan/? = 2 (~ 780 GeV); this also implies that mH+ and uih are 
kept (approximately) constant. However, the b and f masses are still allowed to vary with 
tan/5 as predicted by SUGRA, and the Yukawa contributions to the neutralino couplings 
are included; we see that this suffices to reduce Qh'^ by approximately a factor of 3 at the 
largest allowed value of tan/5. The reduction of rrif^ is the dominant effect here, since for 
the given choice of m and M the b squarks are almost twice as heavy as the f sleptons, and 
have smaller hypercharge. Because for the given choice of parameters one has -C m^, 
the annihilation cross section is essentially proportional to J2i '"^r here, where i is 
a generation index; at small tan/?, all three generations contribute almost equally, but 
for large tan/5 the sum is dominated by the contribution from the 3rd generation. The 
reduction of Qh"^ by a factor of 3 then corresponds to a reduction of by only a factor 
of 1.63. An even larger reduction of the f mass is not possible here since a further 
increase of tan/5 would lead to mp < 0; as shown in ref. this constraint implies 
|tan/5| < rritiQo) / mb^Qo) . Finally, the short dashed curve has been obtained by ignoring 
all effects from the b and r Yukawa couplings; it is a few % above the other curves even 
at small tan/3 since the effects of sfermion mixing are not entirely negligible even here. 
We see that in this case, which approximates the usual analyses based on a global SUSY 
model with independent sparticle masses at the weak scale, VLh? does indeed only depend 
very little on tan/5. 

Using the insight gained from figs. 1-4 it is now quite straightforward to interpret the 
contour plots of figs. 5 and 6. Each of these figures is for fixed values of m and rrit and 
for a given choice of the signs of M and tan/5. In figs. 5 a-d we choose m = 250 GeV, rrit 
= 140 GeV, and explore the plane of M and tan/5 for all 4 combinations of signs. Solid 
and long dashed lines are contours of constant Qh"^ = 1 and 0.25, respectively. The short 
dashed curves in fig. 5a are lines of constant Qh^ = 0.025; since these contours cluster 
very narrowly around the Z, h and P poles, in the other figures we have merely indicated 
the position of these poles with the short dashed lines. 

Finally, the region outside the dotted curves is excluded by the experimental and 
theoretical constraints discussed at the end of sec. 2c. In the region of small \M\ the 
most important constraints are the LEP search limits (29a)-(29d) as well as the gluino 
mass bound (29f); for small |tan/5| the LEP limit on the production of the light scalar 
Higgs boson also plays a role. The lower bound on | tan/5 1 is determined by the Higgs 
search bound as well as the requirement that the Higgs potential should be bounded from 
below even at the GUT scale, as discussed in sec. 2c below eq. (|30|) . The region of large 
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I tan/3 1 and small and moderate values of \M\ is limited by the LEP bounds on associate 
hP production, which is practially equivalent to requiring mp > here, since the overall 
mass scale (m) is chosen quite high in these figures. Finally, in the region of large \M\ 
(or, more accurately, large \M/m\) the parameter space is limited by the requirement of 
a neutral LSP, eq.(29f); this bound is especially important for large |tan/?|, since nif^ is 
smaller there, as discussed in connection with fig. 4. 

The gross features are the same in all 4 figures. At large |M|, VtK^ exceeds 1, as already 
shown in figs. 1; this also happens at small |tan/3| and small |M|, below the Z and h poles. 
Finally, there is a third region where the x I'^lic density is unacceptably large, covering 
the region between the Z and h poles and the VV, Zh and hh thresholds at small and 
moderate values of |tan/3|. The extension of these last two excluded regions does depend 
on the signs of M and tan/5, however. We have already seen that the contribution from 
h exchange and the annihilation cross section into the hh final state are much smaller if 
Ml ■ /i ■ tan/5 < 0, fig. 5b, c; in this case the line fi/i^ = 0.25 can even cross the h pole. 
Moreover, for \M\ < 300 GeV negative values of tan/5 (figs. 5b, d) require rather large, 
positive values of A; if in addition M > 0, X2 of eq.(p8|) and hence become large, 
leading to small higgsino components of the LSP and thus large values of Qh'^ in fig. 
5b. In fig. 5d we have M < 0, however, which results in a partial cancellation in X2 
and much smaller values of together with sizable contributions from h exchange and 
hh production this explains the smallness of the cosmologically excluded region for this 
choice of signs. On the other hand, small or moderately positive values of tan/5 imply 
A ~ 0.5 in the region \M\ < m, reducing the impact of this parameter on X2 and n; 
furthermore, choosing M > now results in larger contributions from the light Higgs 
boson. The differences between figs. 5a and 5c in the experimentally allowed region are 
therefore smaller than those between figs. 5b and 5d. 

Because of the cancellation in the neutralino mass matrix discussed in connection with 
figs. 1, LEP constraints from neutralino searches lead to a more stringent limit on \M\ 
if Ml ■ /i ■ tan/3 > 0; this explains why the experimentally allowed, but cosmologically 
excluded region of small \M\ is larger in figs. 5b, c than in 5a,d. Finally, as shown in 
ref.[^, one can only achieve |tan/5| ^ 1 for sizable values of if A > 0, and usually 
A > m; choosing M < then reduces the effect of the Yukawa coupling on the running of 
scalar masses, as exemplified by X2, eg. (^81). Somewhat paradoxically, this reduces mp, 
because in the relevant limit tan^/5 S> 1 eqs. and imply mp oc X2. A smaller X2 
also implies smaller fi and hence less stau mixing, which again increases m^^. Altogether 
we thus see that choosing M < leads to smaller values of mp and larger m^^ in the 
region of large |tan/5|. The region of parameter space to the right of the P pole that is 
allowed by the requirement m^ < m^^ is therefore somewhat larger for M < 0, figs. 5b, d, 
than for M > 0, 5a, c. One even finds another small region with Qh'^ > 0.25 at very large 
|tan/5| and M ~ -800 GeV. 

In figs. 6 we have chosen M and tan/? to be positive, and study the effects of varying 
rrit (figs. 6a,b) or m (6c, d). We have already seen that increasing m^ increases \fi\, and 
thus also Qh'^ if x is bino-like. Indeed we find larger cosmologically excluded regions in 
fig. 6a than in fig. 5a. Moreover, the fraction of the plane with M < 200 GeV where 
Qh"^ > 0.25 is now much larger than before. This is partly due to the increase of the mass 
of the light Higgs boson caused by the increase of A22 in eq. (p6D . For tan/5 > 5 the h 



23 



and Z poles are now sufficiently far apart to allow for a new region with cosmologically 
interesting DM density in between these poles. Larger values of also imply |^ more 
f mixing, which reduces m^^; at the same time increasing implies larger values of mp, 
see eg. (|25|) . The neutral LSP constraint (29f) therefore does no longer allow to choose 
|tan/3| so large that mp ^ 2m^, except for a small stretch at M ~ 130 GeV; the effect of 
the reduction of mp and m^^ at large tan/3 is nevertheless still quite pronounced in fig. 
6a. Finally we mention that the requirement that the top Yukawa coupling remains finite 
up to scale Mx implies |tan/?| > 2 for m^ = 170 GeV. 

Reducing m^ from 140 to 110 GeV, fig. 6b, has obviously the opposite effect as increas- 
ing it to 170 GeV. is now so small that the whole region of small and moderate \M\ 
is cosmologically safe; however, this also implies larger couplings of the lighter neutralino 
states to the Z boson, so that the LEP search limits rule out a much larger part of the 
plane than in fig. 5a. (Recall that the signs of the parameters are such that cancellations 
occur in the determinant of the neutralino mass matrix.) We already saw in fig. 2 that 
the hh contribution depends very sensitively on M, /z and tan/3; in the given case it is 
large enough to create a small "island" with fi/i^ < 0.25 at tan/? ~ 2, M ~ 200 GeV. 
Moreover, for given values of M and tan/?, ruf^ is now larger and mp is smaller than in fig. 
5a, so that a sizable region of parameter space to the right of the P pole is again allowed, 
including a substantial region where fi/i^ > 0.25. The relatively light pseudoscalar here 
even affects the cosmologically excluded region at very large M, leading to a much steeper 
slope of the uppermost solid line than in fig. 5a. Finally, the reduction of also implies 
that the requirement (^) of a bounded Higgs potential at scale Mx now excludes a large 
region of parameter space at large \M\ and small tan/3. 

In figs. 6c, d we have again chosen vnt = 140 GeV, but have varied m compared 
to the value of fig. 5a. A larger m means larger sfermion masses and hence smaller 
contributions from sfermion exchange diagrams, which are dominant at large |M|, as we 
have seen above. Indeed, we see that the uppermost contour with Qh"^ = 1, which (at least 
for small |tan/3|) occured essentially in the same place in the previous 6 figures, depends 
very strongly on m: For m = 400 GeV the cosmologically safe region above the VV, Zh 
and hh thresholds has disappeared completely; on the other hand, for m = 125 GeV, fig. 
6d, all experimentally allowed combinations of M and tan/3 give Qh'^ < 1, so that the 
requirement that the relic LSP density should not overdose the universe doesn't constrain 
the parameter space any further. Recall that an increase of m also implies an increase of 

via eq.(p7D; this explains the differences in the region of small M, below and around 
the h pole, between figs. 5a, 6c and 6d, even though sfermion exchange contributions 
are essentially negligible here. Finally, for m = 125 GeV we observe a region with very 
small relic density, Qh"^ < 0.25, in approximately the same place that leads to flh^ > 1 for 
m = 250 GeV; this once again demonstrates the strong dependence of the contributions 
from gauge and Higgs final states on the parameters of the model. 

This concludes our discussion of samples of the parameter space of the model. We now 
attempt to derive bounds on sparticle masses or model parameters from computations of 
the DM density. 
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4. Bounds 



We have seen in the previous section that very heavy LSPs tend to have small annihilation 
cross sections. This is not surprising, since unitarity requires the cross section to fall off at 
least like l/w^^ as —>■ oo, for fixed values of the other parameters. Indeed, we see from 
figs. 5 that for given m and rrit the requirement Qh"^ < 1 imposes an upper bound on \M\. 
Moreover, in the important case of a bino-like LSP away from all {Z and Higgs) poles, 
the resulting bound depends essentially only on m and m^, as shown in figs. 1, 5 and 6. 
The SUGRA relations for sfermion masses ( p2|) imply that the sfermions with the largest 
hypercharge, the superpartners Ir of the right-handed leptons, also have the smallest 
masses; they will therefore dominate the annihilation cross section unless -C m^, as 



already pointed out in the previous section. From eqs.(P^) and ( PT| ) and the numerical 
result of figs. 1 that Qh"^ ~ 1 for m = 300 GeV, Tn^= 180 GeV if x is bino-like, one then 
derives the approximate bound 
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(34) 



this bound is only valid for a bino-like LSP away from poles. We already mentioned in 



the previous section that for fixed m the l.h.s. 
Plugging this into the bound (|51D gives 
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has a minimum at 



m < 300 GeV 



0.6m. 



(35) 



for any value of m^; this can also be read off figs. 1. Similarly, for given the l.h.s. of 
(|3^ ) is minimized by choosing m as small as allowed by the constraint m^-^ > m^, i.e. for 



+ 0.83mE = m^,. This immediately gives the bounds 



< 350 GeV; 
\M\ < 825 GeV, 



(36a) 
(36b) 



for any m. The bound (36a) is considerably stronger than the bound of 550 GeV given 
in ref. [0, because in that paper all sfermions were allowed to have mass 



m 



X' 

which is not possible in minimal SUGRA. In particular, a light stop greatly enhances the 
annihilation into ti pairs, which also makes a sizable contribution to the s— wave, 0{v^) 
cross section. 

We emphasize again that the bounds (0) - (^) only hold for a bino-like LSP away 
from poles. In particular, the bound (|35D on m can be violated even for small |tan/?| 
if m^ is close to Mz/2 or m/j/2, see fig. 6c; in this case no useful bound on m can be 
given. Moreover, figs. 5 and 6 show that for a given combination of m, mt and tan/?, 
the upper bound on \M\ is weakest in the region close to the P pole, if m > 150 GeV. 
More precisely, \M\ reaches its maximal allowed value at the point where the contour for 
Qh'^=l (which has a positive slope with increasing |tan/5|) meets the line mf^=mj(, (which 
has a negative slope). The location of this crossing point obviously depends on how large 
TTT'Ti, iTT'P and the Pxx coupling gp^x ai^e for a given choice of M and tan/3. We have 
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already mentioned in the previous section that both \fi\ and mp grow with increasing X2, 
eg. (p8D ■ Moreover, choices of parameters that maximize X2 also maximize the reduction 
of the diagonal elements of the f mass matrix due to the effect of the r Yukawa coupling 
on the RGE. Finally, increasing further reduces m^^ by increasing tl — f/j mixing, 
and reduces gpxx by reducing the higgsino component of x, see eq.(6a). Combinations of 
parameters that increase X2 therefore push the P pole to larger values of |tan/3|, and also 
make it narrower and shallower; at the same time they strengthen the upper bound on 
I tan/3 1 that results from the requirement m^^ > m^. The combined effect of these changes 
is to move the crossing point of the lines mf^ = and Qh'^= 1 to smaller \M\ and smaller 
|tan/3|. 

We already discussed in connection with figs. 5 that in the region of large |tan/?| one 
has A > 0, so that X2 is smaller for M < than for M > 0. The absolute bound on 
\M\ for given m and rrit is therefore reached if M and tan/? are both negative, see fig. 
5d. The resulting upper bound on \M\ as a function of m is shown in fig. 7 for three 
different choices of m^. For m < 125 GeV the upper bound on \M\ only comes from the 
requirement m^^ > m^; for these small values of m the constraint (0) is satisfied even 
for the extreme case = mr . Here the maximal allowed I Ml occurs at rather small 
values of |tan/?|, where nif^ ~ rrif^ only depends on m and M; for small m the bound 
on \M\ is therefore almost independent of nit. However, once the condition Vth? < 1 
starts to impose nontrivial constraints on the allowed parameter space, the bound on \M\ 
comes from the region of large | tan/3 1 and does depend quite strongly on rrit. For light top 
quark and thus small X2 (short dashed curve), the DM constraint reduces |M|max only 
marginally from the value of 5.7m that follows from the simple requirement mj^ > m^- 
On the other hand, if the top quark is very heavy (long dashed curve), a relatively large 
ratio \m/M\ is needed to get sufficiently close to the P pole without reducing m^^ below 
m^. Finally, |M|jnax also depends quite sensitively on the chosen combination of signs; 
for instance, if M and tan/? are both positive one finds M <2150 (1210, 620) GeV for m 
= 500 GeV and rrit = 110 (140, 170) GeV, respectively. 

The rise of the curves in fig. 7 cannot persist indefinetely. Right on the pole, for 
mp/2, the annihilation cross section is proportional to (^p /Fp. The decay width 



m 



Fp is proportional to mp, and thus also to m^; moreover, the PxX coupling decreases 
oc 1/m^ if X is bino-like, see eq.(6a). As a result, the annihilation cross section for a 
bino-like LSP on the P pole decreases oc 1/m^. We estimate that VlK^ will be larger 
than 1 even right on the pole if > 3.5 (1.7) TeV for m^ = 140 (170) GeV. If the top 
quark is lighter, one can arrange = mp/2 even for a mixed LSP, i.e. if |Mi| ^ In 
this case gp^x is not suppressed by small mixing angles, which increases the cross section 
by a factor {mx/MzY compared to the case of a bino-like LSP. The relic density of a 
mixed LSP with m^ = mp/2 will therefore only exceed 1 if m^ > (9(100) TeV or so! One 
can hardly speak of "weak scale" supersymmetry if the lightest superparticle is thousand 
times heavier than the weak gauge bosons; in particular, the model can no longer provide 
a solution of the naturalness problem. 

It can be argued that finetuning is needed to achieve |tan/?| ^ 1, because this occurs 
only over a narrow range of values of A. Moreover, in minimal supersymmetric SU{b) 
the proton decay width increases ||58| oc tan^/3, since in this model proton decay is mostly 
mediated by the exchange of superheavy higgsinos, whose Yukawa couplings grow oc 
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|tan/3|; large values of |tan/?| are then disfavored (although "accidental" cancellations 
might still lead |^8[ to an acceptable nucleon lifetime even if |tan/3| ^ 1). Fig. 7 therefore 
also includes curves (dotted) where we have required | tan/3 1 < 15. This rather mild 
constraint does not affect the curve for = 1 10 GeV at all, but for a heavier top quark 
it leads to a substantial reduction of |M|max- Since ~ mp/2 is no longer possible in 
this case, the only possibility to achieve large values of \M\ is to choose parameters such 
that the LSP is higgsino-like. 

We have already seen in fig. 3 that the relic density of such a state is quite small. 
By extrapolation of the curves of this figure it is clear that a higgsino-like LSP is cos- 
mologically safe up to 
and 
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2 TeV at least; this has already been shown in refs. 
The LSP will obviously only be higgsino-like if \Mi\ ~ 0.43|M| > |yu| 



but 



this requirement might clash with the constraints imposed by radiative gauge symmetry 
breaking, eqs. (P?]) and (P^. Clearly /i^(Qo) can be minimized by choosing A such that 
X2 is minimal, A ~ — 2.08M. For tan^/3 ^ 1, eqs. (|27|) and ( pH]) then give the following 
lower bound on |/i|: 



/x'(<3o) > 



0.9 



150 GeV 
2 



,150 GeV 



2.7-1.04 1 



nit 



190 GeV. 



0.52 



(37) 



Notice that the coefficient of is positive for nit > 85 GeV; it surpasses (Mi/M)^ ~ 0.18 
for rrit > 95 GeV, i.e. in almost the entire experimentally allowed region. The condition 



for having a higgsino-like LSP, |/i(QG 



< 



I Ml I, is therefore most easily fulfilled if |M| 



(and thus |Mi|) is itself very small. However, the coefficient of in eq.(37) also turns 
positive for rrit > 158 GeV. For such a heavy top quark a higgsino-like LSP can therefore 
nothe realized in minimal SUGRA. For this reason the lower dotted curve in fig. 7, which 
is valid for rrit = 170 GeV, is essentially just given by the bound (|34D in the region where 
DM constraints are relevant and the condition ( P5| ) is fulfilled; for m > 300 GeV this 
curve simply follows the h pole, ~ 0.43|M| ^ nih/2. 

For nit = 140 GeV (upper dotted curve) a higgsino-like LSP is still possible provided 
\M/m\ < 0.5 or so. The actual bound on \M\ is in many cases substantially above 
this value since for moderately large m the total XX annihilation cross section is still 
sufficiently large even if the higgsino component of x is sub-dominant; in particular the 
contribution from the ti final state plays an important role here. However, increasing m 
decreases the sfermion exchange contribution to the total annihilation cross section; this 
has to be compensated by increasing the higgsino content of x, i-e. by reducing |/^/M|, 
which via eq.(37) reduces the bound on \M/m\. As a result the bound on the absolute 
value of \M\ is almost independent of m for 400 GeV < m < 1 TeV; for even larger 
values of m the LSP can be an almost pure higgsino even for M > 0.5 TeV, and one has 
|M|max — 0.6m. Finally, we mention that a further sharpening of the bound on |tan/3|, 
e.g. requiring | tan/5 1 < 5, would suppress the bound on |M| for mt = 140 GeV to a value 
very close to the one for mt = 170 GeV. This is because solutions with /i^ <^ always 
have taii^P 3> 1, unless one has \A\ ^ m; however, eqs.(^) and (pSD show that large 
|y4/m| require to be larger than \Mi\, unless mt is near its present lower bound ||59[ of 
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91 GeV. 

This completes our discussion of possible upper bounds on mass parameters of the 
model that can be derived from the requirement Qh"^ < 1. What about lower bounds? In 
principle the bound (p4|) also implies a lower bound on the mass of a bino-like LSP if m is 



fixed within the region allowed by the limit (|35D . Recall, however, that these bounds are 
only valid if x is a nearly pure bino and is not close to a pole. Both these conditions 
can be violated quite easily especially if is not large. Indeed, figs. 5 and 6 show that 
for most combinations of m, and tan/3 the lower bound on \M\ is determined from 
laboratory search limits alone, the exception being the region of small |tan/5| and m > 200 
GeV. In particular, inclusion of the DM constraint does not strengthen the bound > 



20 GeV that follows ||^ from the combination of the gluino, chargino and neutralino 
search limits. 

On the other hand, we have seen in fig. 6d that the LSP relic density becomes 
quite small if m is small, for all experimentally allowed combinations of the remaining 
parameters. Imposing an upper bound on m therefore leads to an upper bound on Qh"^; 
conversely, requiring the neutralino relic density to be larger than some minimal value 
can give a lower bound on m. We have already seen that for fixed m the annihilation 
cross section of a bino-like LSP becomes small both at very small and at very large m^. 



In the experimentally allowed region > 20 GeV and for m < 140 GeV, Qh^ is 



maximized if \M\ is chosen as large as is allowed by the condition nif^ > rriy,. Since m^^ 
decreases at large |tan/?|, Qh'^ will be maximal at small values of |tan/?| where all three 
SU{2) singlet sleptons are essentially mass degenerate. On the other hand, mf^= 
allows ^ m^, which means that all other sfermions will be too heavy to contribute 
significantly to the annihilation cross section, so that eq.(PTD applies. Normalizing the 
cross section from numerical results of fig. 1 as before, we find 

where the constant term comes from the term contribution to mj , which can be 
significant for small values of m; we have checked numerically that the true bound deviates 
from ( pSD by only 10% or so. For m > 140 GeV the maximum of flh"^ for fixed m is reached 
if is at its experimental lower bound; however, in this region eq.(P^) already allows 
Qh'^ > 1 anyway. 

As discussed above, the bound ( p8|) is saturated if m^= m^-^, which implies m^= 2.42 
m OT M = 5.7m; all sparticle masses would then be substantially larger than m. We have 
therefore also studied the question how an upper bound on a physical sparticle mass affects 
the upper bound on Qh"^. We find that fixing the mass of the gluino, of the lighter stop, 
or of the lightest neutralino or chargino state does not induce a significant upper bound 
on Qh'^. On the other hand, we have seen repeatedly that SU{2) singlet sleptons have an 
important effect on the LSP relic density. Since bino-like LSPs annihilate predominantly 
via Ir exchange, a light Ir implies a small neutralino relic density; this has already been 
observed by Roszkowski ||33[| . 

If X is a pure bino, the annihilation cross section (^) is minimized and Qh"^ is maxi- 



mized for given mr if rriy is as small as experimentally allowed EHf, = 20 GeV. For 



small and correspondingly small \M\ the mass splitting between sfermions need not 
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be large; the contributions from all // final states will then have to be included in eg . , 
properly weighted with the fourth power of the hypercharge of the exchanged sfermion, 
and with the sfermion masses given by eq. (p2D . The resulting prediction for the maximum 
of Qh^ as a function of m^-^ is shown by the dotted curve in fig. 8. The other curves in 
this figure are results from numerical scans of the entire allowed parameter space, using 
the program MINUIT of the CERN program library. 

We see that the extended version of eq.(|3T|) does describe the overall trend of the 
numerical results; however, the deviation from the full numerical bounds can be as large 
as a factor of 2 here. In particular, for small values of rrij^ Vth? can be substantially larger 
than one would expect for a pure bino. This occurs if x is a light, strongly mixed state. 
Such a light state exists if \M\ and are both 0{Mz) and Mi ■ /i ■ tan/3 > 0. In this 
case the SU{2) and f/(l)y components of x can have opposite signs (unlike the familiar 
case of a photino); moreover, the squared bino component of this state only amounts 
to typically 50%. As a result, its couplings to charged sleptons are strongly suppressed 
On the other hand, such a state always has sizable higgsino components and thus 



couples to gauge and Higgs bosons. For mt > 130 GeV and present experimental bounds 
on sparticle masses the suppression of the slepton masses is the more important effect if 
mj^ < 150 GeV; since the existence of such a light mixed LSP also implies [|3| a rather 
light chargino the actual numerical value of the upper bound on Qh"^ is in this region 
largely determined by the LEP chargino and neutralino search limits (29a) and (29c). 
For mj^ > 150 GeV the enhancement of the Zxx hxx couplings over-compensates 
the suppression of the couplings to sleptons; in this case choosing x to be bino-like does 
indeed maximize fi/i^, and the full numerical result comes out quite close to the simple 
approximation based on eq. (pi|) . 

The curve for = 110 GeV looks quite different from those for larger values of m^. 
One reason is that for m^^ < 75 GeV the Higgs search limits now force one to choose 
I tan/? I substantially larger than 1; this increases the Zxx coupling, which vanishes for 
|tan/?| = 1. In this region VLh"^ is therefore again maximized by choosing x to be bino- 
or photino-like. For larger values of frti^-, the constraints imposed by the Higgs search 
limits are less severe, since the Higgs boson masses tend to increase with the overall SUSY 
breaking scale; for 80 GeV < m^-^ < 130 GeV the curve for rrit = 110 GeV is therefore 
close to those for a heavier top quark. However, for even larger values of ^[^, ^h"^ is 
maximized for such small values of \M/m\ that radiative gauge symmetry breaking with 
a light top quark can only be achieved if l^/ml is quite large, which greatly reduces the 
mass of the lighter stop eigenstate; the bound m^^ > 45 GeV then rules out large regions 
of parameter space, causing the curve for rrit = 110 GeV to again fall below those for 
larger values of if m^-^ > 130 GeV. 

In all cases we find that the upper bound on Qh"^ for light is determined not only by 
mj^ itself, but also by the experimental bounds on the masses of the other sparticles and 
Higgs bosons. This dependence is further illustrated by the long-short dashed curve in fig. 
8, which is valid for rrit = 140 GeV after we impose the (hypothetical) bounds ^tx,x+ — 
80 GeV and rUg > 200 GeV. This is meant to approximate the bounds that would emerge 
if LEP200 and the tevatron fail to find evidence for supersymmetry (other than perhaps 
light sleptons). The contraints from searches for Higgs bosons and neutralinos will also 
become more severe in the next few years, but the final bounds depend quite strongly 
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on the energy and luminosity that will be achieved in the LEP upgrade. The increased 
bound on the chargino mass excludes a mixed LSP unless > 40 GeV, which is 

quite close to the Z pole; the upper bound on flh'^ is therefore now always saturated if x 
is a bino- or photino-like state. The increased lower bound on the gluino mass implies 
rri)^ > 35 GeV for such a state, see eq.(^); this enhances the annihilation cross section 
by approximately a factor of 2, compared to the present bound > 20 GeV. Since x 
is now always gaugino-like, the bound that can be derived from eq.(^) reproduces the 
long-short dashed curve to an accuracy of about 20%. 



5. Summary and Conclusions 

In this paper we have computed the relic density Qh"^ of LSPs produced during the Big 
Bang, within the framework of minimal N = 1 Supergravity models with radiative gauge 
symmetry breaking. In sec. 2 we briefly described the calculation of the relic density for 
a given LSP annihilation cross section. We then discussed the various contributions to 
this cross section, including final states composed of gauge and Higgs bosons that become 
accessible for a "heavy" LSP, > M^r. Expressions for all these cross sections are given 
in the Appendix, for a general neutralino eigenstate; to our knowledge such a complete 
list does not exist in the literature. In sec. 2b we discussed the qualitative features of the 
most important contributions to the annihilation cross section for the important special 
cases that the LSP is an almost pure higgsino or bino. In particular, we pointed out that, 
even though the higgsino component of a bino-like LSP vanishes as oo, the matrix 

elements for xx ~^ W ^ind XX ~^ remain finite, where V = W, Z and if is a Higgs 
boson. This is because of the enhanced production of longitudinal gauge bosons; the same 
result can also be derived from the equivalence theorem. 

In sec. 2c we briefly summarized the constraints imposed on the particle spectrum 
by the assumption of minimal supergravity (SUGRA). The relations between the masses 
of the gauginos and those of the superpartners of the light quarks and leptons are by 



now well known, and have been included in several previous analyses ^ ^ p4| . 
However, SUGRA also implies relations between the mass of the top quark and the soft 
breaking parameters on the one hand, and the masses of the Higgs bosons as well as the 
supersymmetric Higgs(ino) mass parameter /i on the other; these relations had previously 
only been included |2^ for the case that the SUSY breaking gaugino mass M is much 
larger than the scalar mass m at the GUT scale. We saw in our numerical examples of sec. 
3 that these relations can have large effects on the relic density. In particular, we found 
that the coupling of a bino-like LSP to gauge and Higgs bosons is usually suppressed 
for a heavy top quark and/or large m, since increasing rrit or m tends to increase 
which in turn reduces the higgsino component of x- This strongly affects Qh"^ both in the 
pole region [rriy^ ~ Mz/2) and in the threshold region (m^ ~ Mz)- Moreover, as already 
pointed out in ref. ||3^, flh'^ is greatly reduced if the ratio tan/? of the vevs of the two 
Higgs bosons is large, since this implies large b and r Yukawa couplings, and thus a light 
pseudoscalar Higgs boson P and reduced masses for the light b and f eigenstates. 

The relations between particle masses and neutralino couplings imposed by SUGRA 
also imply that the region of parameter space where the LSP density can lead to a fiat 
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universe (0.25 < Qh'^ < 1), as predicted by inflationary models, has a very complicated 
shape, since the XX annihilation cross section depends on all input parameters; the same is 
true for the region where LSPs can at least provide the DM halo of galaxies {Qh"^ > 0.025). 
Nevertheless it is clear from the results of sec. 3 that a cosmologically interesting relic 
density is obtained quite naturally, provided that x is gaugino-like and the SUSY breaking 
parameter m is not too small. If x is higgsino-like and > Mw, XX annihilation into 
pairs of gauge bosons is so strong that cosmologically interesting values of Qh"^ only 
occur |19| for > 0.5 TeV. Since in minimal SUGRA the gluino as well as most 
squarks are at least 6 times heavier than the LSP, requiring > 0.5 TeV leads to quite 
severe fine tuning |]35| , |16| and is thus unattractive. The xX annihilation cross section for 
higgsino-like x can be quite small below the WW threshold if the gaugino mass \M\ is 
large; however, using results of ref. we estimated that in this case co-annihilation of 
X with the next-to-lightest neutralino x' will suppress flh"^ to a value well below 0.25, 
since the off-diagonal Zxx' coupling is large, and since co-annihilation can proceed from 
an s— wave initial state. We therefore conclude that the higgsino does not make a very 
attractive DM candidate in minimal SUGRA models. 

Unfortunately we saw in sec. 4 that the complicated shape of the region of parameter 
space that yields Qh"^ < 1 makes it difficult to derive interesting upper bounds on mass 
parameters from the requirement that relic neutralinos do not overdose the universe. In 
particular, m can be almost arbitrarily large, and thus all sfermions and most Higgs bosons 
can be very heavy, if ~ Mz/2 or ~ 'm./i/2, where h denotes the light scalar Higgs 
boson. In this case the gluino, the lighter chargino and the next-to-lightest neutralino 
would still lie in the region accessible to the next generation of accelerators {rrig < 500 
GeV, m^+^y < 150 GeV). However, if |tan/3| is large, mp can be greatly reduced so that 

can be close to mp/2. Since P exchange mediates xX annihilation from an s— wave 
state, Vth? is greatly suppressed in this case, allowing both m and \M\ to be well beyond 
1 TeV without getting in conflict with cosmology. If the top quark is not very heavy, 
mt < 155 GeV, cosmologically safe solutions with very large m and \M\ also exist for 
small |/i|, so that the LSP is higgsino-like. In this case one would still have a "light", 
almost degenerate SU{2) doublet of higgsinos, but all other sparticles would be very 
heavy. Since the mass splitting could be as small as C(l GeV), such light higgsinos would 
probably be very difficult to observe even at e~^e~ colliders [^. This scenario is therefore 
very similar to the case of a very large as far as collider experiments are concerned. We 
have to conclude that the requirement Qh"^ < 1 does not strictly exclude the possibility 
that sparticle masses are well beyond the reach even of the planned supercolliders. 

On the other hand, interesting upper bounds can be obtained for the "generic" case 
of a bino-like LSP away from poles. We found that in this case the total annihilation 
cross section is dominated by the final state produced via Ir exchange (/ = e,fi,T). 
Contributions from II exchange and q exchange are suppressed by the larger masses 
and smaller hypercharges of these sfermions, while the contribution from Z exchange is 
suppressed by the small Zxx coupling. We mentioned earlier that the VV and Zh final 
states also contribute with full gauge strength if the LSP is sufficiently heavy, but the 
relevant hypercharge here is the one of the Higgs bosons, leading to a suppression factor of 
1/16 compared to the Ir exchange contribution. This allowed us to derive the very simple 
and yet quite accurate expression ([3T|) for the total annihilation cross section, leading 
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to the analytically derived bounds (^) and (pq). In particular, we find for this case an 
upper bound rrig < 2 TeV, which is only slightly less stringent than the bound derived 
[p5| , p!6| from the requirement that fine tuning should occur at most at the 10% level; 
the corresponding bound on is considerable more stringent than the one that follows 
under similar assumptions [jl8[ in a more general SUSY model. Our bound m < 300 GeV 
is potentially even more interesting, since this value is below the one derived |3^, 16| from 



fine tuning arguments; together with the bound on rrig it would virtually guarantee that 
at least the sleptons, the light chargino and the next-to-lightest neutralino should be 
observed at a TeV e'^e~ supercollider. However, we have to remind the reader that there 
are several ways to evade these bounds. 

As first pointed out in ref. and also emphasized in refs. [p3|, p^ , the LSP relic density 
can only lead to a flat universe with 0.5 < h < 1 if m is not too small. We quantified this 
in the simple relation (P^, which shows that Qh"^ > 0.25 is only possible for m > 40 GeV. 
(In ref. the more stringent bound m > 100 GeV has been found, but this is only valid 
for < Mw-) Furthermore, requiring Qh'^ > 0.25 also implies a lower bound on the 
mass of SU{2) singlet sleptons. The exact value of this bound depends on the top mass, as 
well as on the bounds on the masses of other sparticles and Higgs bosons, but it is always 
close to 100 GeV, see fig. 8. We thus see that requiring 0.25 < Qh"^ < 1 determines 
mj^ to be within a factor of 2 of 200 GeV, if the LSP is gaugino-like and not near a 
pole. Unfortunately, in this case LEP200 will fail to discover a slepton. This has already 
been pointed out in ref. for a light LSP, but at least within the framework of SUGRA 
models this conclusion also holds > Mw- (Of course, we need < always.) 
On the other hand, the gaugino mass parameter M is only poorly determined by the 
requirement that Vth"^ lies in the cosmologically interesting range, since the annihilation 
cross section is almost independent of M over a wide range, see figs. 1. In particular, the 
gluino, the squarks, the lighter chargino and at least one neutralino could all lurk "just 
around the corner" , but their masses could also lie in the range that can only be covered 
by supercolhders. 

Our overall conclusion is that, while limits on the relic neutralino density allow to 
rule out large regions of parameter space, they do not allow to derive upper bounds on 
sparticle masses which are both interesting for experiments at existing or planned colliders 
and valid for all combinations of the other parameters. On the other hand, results for the 
perhaps most natural case of a gaugino-like LSP do indicate that sparticle masses should 
lie in the range to be covered by planned e^e~ and pp supercolliders. Moreover, if the 
dark matter in our galaxy does indeed consist of neutralinos, one would expect their mass 
to lie within the range of sensitivity pl| of next-generation direct search experiments and 
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bound on the expected signal size is again difficult to derive. Cosmological arguments can 
therefore supplement and lend support to direct SUSY searches at collider experiments. 
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Appendix A: Matrix elements 

Since we are interested in the nonrelativistic hmit of the XX annihilation cross section, 
we employ the partial wave formalism. In this formalism the helicity amplitude for the 

process 

x{h) + x{h)^X{\x) + Y{\y) (A.l) 
h, \x and Ay are the helicities of the corresponding particles) is expanded as follows: 

oo 1 L+S 

T=Y.Y. E (A.2) 

L=OS'=0 J=\L-S\ 

Here the reduced partial wave amplitude A describes annihilation from an initial state 
with definite spin S and orbital angular momentum L, and thus also with definite C and 
P quantum numbers. The spin projectors V depend only on h and /i, and the angular 
dependence is contained in the d functions c^a. a^- \i = h — h and A/ = Ax — Ay are 
the differences of the helicities of the initial and final particles, respectively. Because our 
initial state consists of two identical Majorana fermions, we only need to consider initial 
states with C = 1. Furthermore, since we want to expand the total annihilation cross 
section only up to 0(f ^), only annihilation from s— and p— wave initial states has to be 
included. Altogether we thus find that we need to include only the contributions from the 
^5*0, ^Pq, ^Pi and ^P2 initial states; explicit expressions for the relevant V can be found 



in ref. |^3[. The annihilation cross section is then given by 

a{xx - XY)v = [\A{'So)\' + ^ [\A{^Po)\' + WPi)\' + \A('P2)\'] } • (A.3) 

Here v is the the relative velocity of initial neutralinos and s is the total energy. S is 
a symmetry factor which is 2 when X = Y. The summation over the final helicities is 
implicit in this equaiton. Finally, the kinematical factor /5/ is given by 




f3f = \l- 2 ^ + ■ (A.4) 



In this Appendix we list the helicity amplitudes A(^^~^'^Lj) for all two-body final states 
accessible to xX annihilation in leading order in perturbation theory. We first list some 
couplings which appear in many matrix elements: 

O^^ = -^No4V,2 + No^V,^ (A.5a) 
= -^No3Uj2 + (A.5b) 



1 

V2' 



Qo5 = ^04^,1 + -7^{Noj + A^oi tan 9w)V,2 (A.5c) 
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Qoj = NosUji - ^(iVo2 + iVoi tan ew)Uj2 (A.5d) 



V2' 

2- .3+ 2- 



1 



Ooj = -o^03A^,3 + T^No^N^^ (A.5e) 



Qoj = o [A^03(A^,2 - tan^^H^iV,i) + (0 ^ j)] (A.5f) 



•^o, = 2 [A'o4(A^,2 - tangly iV,i) + (0 ^ j)] (A.5g) 

The expressions for O^, and O"^ are the same as in ref.0; our definitions for 
Q'^, Q'^, Q" and 5"' differ shghtly from those of ref.[^. U and V are the matrices 
that diagonahze the chargino mass matrix A4^, which is given by [^]: 

^ - 1^ Mh.v^cos/3 /i )■ ^^-^^ 

The matrices f/ and \^ can be chosen to be reah 

y COS0+ sin0+ \ 




In the hmit where either IM2I or is much larger than M^, the chargino mass eigenstates 
are almost identical to the current states. The (small) mixing angles then become: 

^- = iS^{M2Cos(3 + fism(3y, (A.8a) 
M2 — 

0+ = yF^^, (M2 sin /? + /i cos /?) . (A.8b) 

The Nij in eqs. (|A.5|) are elements of the matrix that diagonalize the neutralino mass 
matrix, eq.@; as already discussed in sec. 2b, it also becomes almost diagonal if IM2I 



or is very large, see eqs.(|^). Unlike refs.Q and we do not require the chargino 
and neutralino masses to be positive; Nij can therefore also chosen to be real. We have 
checked by explicit calculation that this still preserves all relative signs in our amplitudes, 
provided we keep the sign of the fermion masses everywhere. Finally, the suffix zero in 
eqs.( |A.5| ) refers to the lightest neutralino, x- 

The coupling of the neutral Higgs bosons to neutralinos can be expressed in terms of 
Q" and S" [gg: 

Tpoj = - sin(3Qlj + cos/JSg^-; (A. 9a) 

^loj = - cosaQlj + sina^oj; (A. 9b) 

T20j =sina(5o'j + cosa^gj; (A. 9c) 

(A.9d) 

Here the labels P, 1, and 2 refer to the pseudoscalar, and the heavier and lighter neutral 
scalar; a in eqs.( [A.9D is the mixing angle of the neutral scalar Higgs bosons as defined in 



ref.[42], including leading one-loop radiative corrections. 
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Finally, we define some kinematical quantities and propagators: 



A 



2 rrix + iTLy 



2ml 



Rx = mx/m^, Ry = niY/m^, Rj = mj/m 



(A.lOa) 

(A.lOb) 
(A.lOc) 



Here, rrij is the mass of an exchanged particle, and P/ the v —>■ limit of the inverse of 
the corresponding t ot u channel propagator. 

We are now ready to list the partial wave amplitudes of the contributing processes. As 
mentioned above, we only need the nonrelativistic limit of the annihilation cross section, 
i.e. its expansion in powers of v up to and including terms of However, this 

expansion breaks down |^ in the vicinity of poles and thresholds, since there the "higher 
order" terms can actually diverge. In our expressions below we therefore only include those 
terms of order that result from the expansion of well-behaved functions of v. Specifically, 
we include terms that result from the expansion of t— and m— channel propagators, as well 
as terms that result from the calculation of spinors or bosonic wave functions. On the 
other hand, we do not expand the kinematical function (3f, nor s— channel propagators. 
Of course, far away from the threshold or pole, these terms can be expanded safely. This 
can be incorporated into our matrix elements by the following substitutionsiB 



Pj^Pj{v = 0) + 



iPfiv = 0) 

„,2 \ 



X 



my) 



A-Rj 4-Rj 



A-Rj 



(A.lla) 
(A.llb) 



Clearly, these substitutions should only be used in the C(f °) terms of the ^5*0 amplitudes. 
Because of the suppression factor S/xi? in the expression for VLh?, as well as the small- 
ness of the the numerical factors in eq.(A.lla), the numerical effect of these substitutions 
is only sizable in cases where the expansion itself can no longer be trusted; in this case 
only the much more complicated methods described in ref.[^] give reliable results. Since 
a more careful treatment of poles and thresholds will not change our conclusions, we do 
not pursue this avenue here. It should be noted that the 0{v) terms that result from 
the expansion of t— and m— channel propagators can, e.g., change the annihilation cross 
section into // final states by a factor of two; fortunately the 0{v) terms are always 
regular. 

We list contributions with different final state helicities separately; all these contribu- 



tions have to be added incoherently, as shown in eq. (|A.3| ) 



1) xx^w-{\)w+{X) 

This final state receives contributions from the exchange of the two chargino eigenstates 
(labelled by j) in the t— or u— channel, as well as from the exchange of the two neutral 



*Far away from the pole, the propagator can be taken to be real. 
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scalar Higgs bosons (labelled by i) as well as the Z boson in the s— channel. In the 
following expressions, summation over subscripts (i, j, . . .) that appear more than once 
is always understood. 



A^So) : A/ = 0, A = ±1 



2R: 



1 _ A ^ 
4 ~ m ^ 3^ 



(A.12a) 



^(3Po) : A/ = 0,A = 



w 



1 - 



3R 



, qR'2 



3R 



3R 



(4-i?|,^+iG'^,jit:M^ 

Xf = 0,\ = ±l 
Ai'Pi): \Xf\ = l, 



20gOgi?;\^ VQvg'RwF, 



4 - + tGn^ 



(A.12b) 



(A.12c) 



w 



+ 2vPfg 



,2 "-^Oj 



qR'' 



RiyPj 



1 - 



2 - 



_20gOg^1 
J 



-Rw(4 — R\) 



= 0, A = 



RwPj 



00 



A/ - U, A - ±1 (Oo^- - Ooj j ^_^2 

A{^P,):\X,\^2 -^/(OgVog.^) 
|A/| = 1 



(A.12d) 

(A.12e) 
(A.12f) 
(A.12g) 
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2vg^ 
R 



w 



""^^ (O^.; + off) + + A,/5/<^ } (A.12h) 



P2 



Pf 



= 0,A = ±1 

|Ay| = 0,A = 

-.2 ( nL"^ 



VSR^ 



ij2 



W 



p. 



p2 

J 



R2 



(A.12i) 



(A.12j) 



Here, the i? are the rescaled masses of exchanged or final state particles, see eq.A.lOc) 
and we have introduced the rescaled widths 



G 



(A.13) 



The definition of the Pj is as in eq.(A.lOb), and the Fi are given by |52 

Fi = cos(/3 - a)Tioo, F2 = sin(/3 - a)T2oo 



(A.14) 



2) XX ^(A)^(A) 

Here the contributing diagrams are very similar to those of the W^W~ final state, except 
that the Z exchange contribution is absent. Moreover, the subscript j now labels the 4 
neutralino eigenstates. The nonvanishing partial wave amplitudes can be written as: 



A{^Sq) : A/ = A = ±1 



f 



+ 



Ai^Pn 



P^ 

A/ 



3P, 



+ 



R^ 



2m 



3Pf 



0,A = 
2 



'Oj 



3R- 



+ 2VQvglO''f 



Rj 




2 " 




Pj 


1 - 


3P,J 


1 


1 - 


^PA 


Rj 
+ — 
P 



{4-Rl^+iGH,)Rw 

Ay = 0,A = ±1 
2v/6.,|0;f U 



3R 



P^ 



2m 



3P, 



y/6vglRwFi 



P 



4 - Rjr + iC 



(A.15a) 



(A.15b) 



(A.15c) 



37 



A('Pi): \Xf\ = h 



oj 





1 ■ 








^1} 



A{^P,): \Xf\=2 - 



|A/| = 1 



R? 



\Xj\=0,\ = 



v/3i?r°^ IP 



+ 



2R.f3 



P2 



In eqs.( |A.15D we have used the usual notation 



9z 



cos9 



w 



(A.15d) 
(A.15e) 
(A.15f) 
(A.15g) 
(A.15h) 



(A.16) 



3) XX - Zi\)H, 

This final state receives contributions from the exchange of the four neutralinos (labelled 
by j) in the t— or m— channel, as well as from the s— channel exchange of the Z boson as 
well as of the pseudoscalar Higgs bosons P. The result is: 



A{^So) : A = 

Rz 



99z 



2Jj{Rj - 1) 



+ 



O 



4L 



RzcosO 



w 



4-RI + iGi 



— v 



Rz~'''T, \ 2 



99z^^^{R,-^) 



2{Rj - 1) A{Rj - 1) - 



+ 



+ (2- 



3P 



2 V ^iJf 

V -^^99z 



3L 



O 



R, 



A-Rl + iGp ARzcosO 



w 



A{^Pi) :A = ±1 



^^99z-^ 
j 



■ 2 ^4 
-^1 T 

1 4 



JjRn 9 ^ Rz 



A = 



27^(1 + 17)^?^?^ 2(1 + P,)-f 



R 



Pj cose w {4: - R 



-O 



A{'P2):X = ±1 -AXvp/ggz^ 



(A.17a) 

(A.17b) 

(A. 17c) 
(A.17d) 
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In addition to A, which has aheady been defined in eq.(A.lOa), we have introduced the 
following quantities: 



Tpmp 



mi 



(A.18a) 







iOL) J 





2 


J, 


= - OlfTaOj for Ha] 


(A.18c) 


L 


= ^ sin(Q; - (3)Tpoj for Hi] 


(A.18d) 




^ cos(a - l3)TpQj for H2] 


(A.lSe) 





=Olf sin{a - 13) ioi Hi] 


(A.18f) 




Olf cos(a - 13) for H2. 


(A.lSg) 



4) XX - Z{\)P 

The contributing diagrams are very similar to those of the ZH final state, except that 
the single P exchange diagram has to be replaced by Hi, H2 exchange diagrams; notice 
that these scalar Higgs bosons cannot be on-shell, so that we need not introduce complex 
propagators here. 



A(^P„) :A = 



99z 



- + i?? \ 

3^36 



3 



Li 



A{^Pi):X = ±l -AvPfXfggz 



-Ri + — 
^ 2 



3 



A{^P2):X = ±1 AvBfggz 



-R 



A = - 
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73 



l + i?,-A2 + -(i?,-l) 



In eqs.( |A.19D we have defined: 



(A.19a) 
(A.lQb) 
(A.19c) 
(A.lQd) 



Li =^ sin(a - /5)rioo, 
^3 =OofTpoj 



^2 = cos(a - /?)T2oo 



(A.20a) 

(A.20b) 
(A.20c) 
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5) XX^W-{X)H+ 

In this case one has contributions from the exchange of the two charginos, which are again 
labelled by the suffix j, as well as from the exchange of all three neutral Higgs bosons. 
Notice also that there are equal contributions from W^H~ and W~H^ production, lead- 
ing to an overall factor of 2 in the final cross section in this case. The contributing partial 
wave amplitudes are: 



A(iSo) :A = 



D__. D ^ D__.D 1 O p. 3p2 Q c 



3R 



2 - 



RwPj 
1 



^ ' RwPj 1 3P2 



4^/_^ 1 



R\Y 4 — Rt 



1 + 



(A.21a) 



Ai'Po) :A = 



w 



1 + 



2 „2 2^2 ^ 

3 3 6 



Rw 

A{^Pi) :A = ±1 



(A.21b) 



- 4f 
4f/?/A 



p2 



+ 



3 ^ J 
2 



2 



A = 
A(3P2) :A = ±1 



RwPj 



{J- + J';) + {.]■- J-)Rj 



P2 



/j;+^9M(/<-;+<)y-«;(a-:;.-a-;; 



(A.21c) 
(A.21d) 

(A.21e) 



A = 

SvPjg^ 



\/?)RwPj \ 
In eqs( |A.21| ) we have introduced: 



-)(^;-^;)} (A.2if) 



=ij{Rw - Rh+) 



(A.22a) 
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- j(05 


-Oo^,)(g;fsin/3 + Q;5 cos/3) 


(A.22b) 




- i(05 


+ Oo^,)(g;fsin/?-Q;5 cos/3) 


(A.22c) 


jy 




- <^oj)(voj sin/3 - Qoj cos/3) 


[A.zza) 




- i(^o1 


+ C'o^,)(QoiSin/? + Q;5cos/?) 


(A.22e) 
(A.22f) 



6) XX ^ ^a^^f, or PP 

Both the HaHh (a, 6 = 1,2) and the PP final state can be produced by the exchange of 
one of the four neutrahnos (labelled by j) in the t— or m— channel, as well as by scalar 
Higgs exchange in the s— channel. The contributing amplitudes for the HaHf, final state 
can be written as: 

A (3Po) -.VQvg^ \9abiT^oo-. „f ^, - 2TaOjnoj^^^ + ^T,o,TM,%lA.23a) 

[ 4 - + iLxH, Pj J> P, J 

A (3P2) : - l=vp/g^TaOjTboj^ (A.23b) 

The corresponding PP amplitudes can simply be obtained by replacing a, b by PP and 
Pj by — Pj. The rescaled width Gj^. = V HiTnUi/Tril^ as before, and the trilinear Higgs 



couplings are |^ 



3 

(7111 = — cos2a; cos(a + /?) (A. 24a) 

2cos6'i4/ 

3 

fi'222 = -t; 7— cos2q; sin(a + /?) (A. 24b) 

2cos6'vy 

9112 = 9121 = — [2 sin2a cos(a + /?) + sin(a + /3) cos2a] (A. 24c) 

2cost'vF 

fi'122 = ^7221 = n — [2 sin2Q; sinfa + /3) — cos(a + /3) cos2a] (A.24d) 

2cos6'v^ 

cos(a + 0) / . X 

^ppi = 9iPP = ^ n COS2/3 A.24e 
2cos6'vi/ 

sin(a + /3) / . ^ 

9PP2 = 92PP = 7^^^ COS2/3 (A.24f ) 

Recall that the suffix 1 (2) refers to the heavier (lighter) Higgs scalar. 
7) XX ^ H,P 

This process receives contributions from diagrams where one of the four neutralinos (la- 
belled by j) is exchanged in the t— or u— channel, as well as from s— channel exchange of 
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Z OT P bosons. The result is: 
A {'So): 

2V2g'^gappTpoo 



Rz 
A-Rl 



1 + 



V2g' 



rip — It, 



AV2g^TaOjTpoj^ i^l + V 
V2{Rl-Rl)g^Tao,TpoAl + 



9aPZ- p2 



2P, + 3Pf, 



1 



4 - i?| + zGz 



2 TaOiTpoi 



(A.25a) 
(A.25b) 



Here, = M^F^/m^ as usual, and we have introduced the combinations of couplings 



sin(a-/5) „i 

gipz = — z — — — u, 



00 



g2Pz 



2 cos'^Ow 
cos{a- (3) „^ 



2 cos2^ 



w 



00 



(A.26a) 
(A.26b) 



8) XX 



This final state gets contributions from exchange of the two chargino states, labelled by 
j, as well as from the exhange of the Z boson and the two scalar Higgs bosons, labelled 
by i. The nonvanishing amplitudes are: 



A('Po): -V6vg'{{Q'of + Q'o'} 



1 2/5j 



P, 3Pf 



'L ' P^ 
^3 



+ VQvg'^Rw 



A (3P0 : 2v(3fg\Q',)^ - Q,f)- + Av(3f 



'f?2^ 1 . . ^ COs29w 2 Oqq 



p. 



3 

2. 1 



cos2^w^ A-R 



(A.27a) 

(A.27b) 
(A.27c) 



In eqs.( |A.27D we have introduced the quantities 

= —^100 

9+-2 = — 7^200 



cos(a + /3) 
cosfp - a) — cos2/y 

2 COS^'fe'iy 



sin /? -a)+ ^ ^ COS2/3 
2 cos^fc'vF 



(A.28a) 
(A.28b) 
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9) XX ^ fa{h)faCh) 

Here we use the suffices a, b to label the up and down components of weak isodoublets, 
with T^^a = ±1/2. This allows us to use the same notation for quarks and leptons. Of 
course, in the final cross section a color factor of 3 has to be included in case of quarks. 

This final state can be produced by the exchange of the two sfermion eigenstates in 
the t— and it— channel, as well as by s— channel exchange of the Z boson or of one of the 
three neutral Higgs bosons. Since the final state particles can be massless, we have to 
include finite widths for all s— channel propagators. We remind the reader that we do not 
expand s— channel propagators with respect to v, while we do expand t— and m— channel 
propagators. The contributing partial wave amplitudes can be written as: 



A{'So) :A/ = 



2Pi ' 3Pf 

1 1 



4 2Pi 6Pi 3P{ 



+ (-1) 



h+ 



1 2V2g'^ 



cos2 ew~'"''''"'R 



^00 ^3a^ 
Z 



+ 4V2(-l)^+25/ip,Tpoo 



4-RI + iG 




A i'Po) :Xf = 



- 2\/6v/3fg 



Pi' 



hiaTiOO 



'4-Rjj^+iG 



Hi 



A (3Pi) :A; = 



vR 



Pi 



HKo' - v^^o') + (Xo - Ko, Ko - Ko, Pi - P2) 



4 - P| + iGz 



(A.29a) 



(A.29b) 



(A.29c) 



V2v xMx'J + W'J 



Pi Pf 



1 1 



+ 2^;A;(XoW^:o)^ 
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+ 2V2vglOl^{+XfT,Jf 



T3,a - 2e/^ sin^ 9w 



} 



4-RI + iGz 



(A.29d) 



A i^P,) -.Xf = 



(A.29e) 



(A.29f) 



In eqs.( |A.29| ), we have defined Rf = rrif^/m^ as in eq.(A.lOc). These expressions fully 
include mixing between SU{2) doublet and singlet sfermions. The sfermion mass eigen- 
states, labelled by 1 and 2 in eqs. (|A.29|) , are defined by 



/i = fLCOsOf + fRsmOf, 
h = -h sin6'/ + Jr cosOf] 



(A.30a) 
(A.30b) 



the corresponding rescaled masses determine the inverse propagators Pi and P2 as shown 
in eq.(A.lOb). Sfermion mixing also affects their couplings to neutralinos: 



X'^Q = Xao cos9f + Zao sm9f, W'^Q = Zao cosOf + Yao sinOf, 
= -^ao sinOf + Zao cos9f, F^'o = -Zao sinOf + Yao cos9f. 

The couplings of the unmixed sfermions in eq.(A.31) can be found in ref . |[42||i: 



Xao = -V2g [T3aNo2 - tanOwiTsa - e/jA^oi] 
Yao = V2gtsin9wef^Noi 

No3 



ZuO 



grriu 



v/2sin/3Miy 
V2 cos/3 Mw 



04 



Finally, the couplings between Higgs bosons and SM fermions are fi^ : 



(A.32a) 
(A.32b) 

(A.32c) 



grriu cot/3 



2M, 



girid tan/5 



w 



2M, 



w 



(A.33a) 



^Notice that the last term in the relevant eq.(5.5) of that reference has a wrong sign; H. Haber, private 
communication. 
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hiu 

h2u 



gm^ sma 
gniu cosa 



h 



2d 



grrid cosa 
~2Mwcos(3 
grud sina 
2Mw cos/3 



(A.33b) 
(A.33c) 



Appendix B: Application of the Equivalence Theo- 
rem 

As already discussed in sec. 2b, the high-energy hmit of the amphtudes for the production 
of longitudinal gauge bosons can be understood from the equivalence theorem [^. This 
also provides a useful check of these amplitudes. As an example, we discuss in this 
Appendix the production of two longitudinal Z bosons. We see from eqs. ( |A.15| ) that only 
two amplitudes contribute to the A = A = final state; in the limit \m^\ ^ Mz they 
become 



RIP, 



I- 



■J 



2y/6vg2Fi 



^vgl 



v/3P|p/°^- P,, 



(A.34a) 
(A.34b) 



where we have already made use of the fact that O'^^ oc so that only contributions 
with an explicit factor 1/-R| survive in the high-energy limit. More specifically, we see 
from eqs.(^ that only the contributions from the higgsino-like neutralinos (j = 3,4) will 
survive in this limit: 



R?. — —Ra 



Ml 



P. = P. = 1 + 



/i 



^03 



+ 0, 



"l2 
04 



Mr 

1 / MzsinOw ^ 
4\m[^j 



(M2 + /i2 + 2Mi/isin2/?) 



•^03 



1 /M^sin6' 



w 



4 V - /i2 



Mf + fx^) sin2/3 + 2Mi/i 



(A.35a) 
(A.35b) 

(A.35c) 

(A.35d) 



In SUGRA, ^ M| almost always implies mp ^ M|. In that case Fi becomes very 
small, so that the contribution from the heavy Higgs scalar can be neglected; moreover, 
one has Rj^^ <^ 1 in this limit, and a = (3 + 7c/2. This gives: 

M^sin^^H^ (M,+/.sin2/3), (A.36) 



200 



COS^VK - Ml 

where we have used eq.(6a). Inserting eqs. ( [A. 351) and ( [A.36[ ) into ( |A.34| ) finally yields: 



A(^Po 



VQvg 



i2 



Ml 



{^M^ [Ml + /i^ + 2Mi/xsin2/? 
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-Mi/i 1 



Ml 



3 Mi + /i2 



(mI + /i^) sin2/3 + 2Mi/i 



Ml 



V6 ,2 

2 ^ V-Mf 



^■'^^^ T72 {Ml +11 sin2/5) ; 



(A.37a) 



M2 



v/3 ^ M2 + /i2 \Ml-ii\ 
2/iMi 



{Mi2 + /i2 + 2Mi/isin2/? 



M/ + sin2/? + 2Mi/x 



(A.37b) 



The last term in eq.(A.37a) comes from Higgs exchange. One recognizes the form of our 
"symbohc" expression (|13D of sec. 2b. 

On the other hand, since the neutral Goldstone boson G is a pseudoscalar, the am- 
plitudes for the production of a GG pair can be directly read off from our results for PP 
production, eqs.( |A.2^ ). Retaining only those terms that remain finite as m'^/M'^ oo, 
one has 



A(^P, 



-^GOj 



P. 



3Pi 



(A.38a) 
(A.38b) 



The couplings Tcoj can be obtained from eq.(A.9a) by the substitution sin/? — > — cos/3, cos/3 
sin/3. We see that again only the contributions from the exchange of the higgsino-like 
neutralinos (j = 3, 4) survive: 



Tgo3^ + Tgo4^ — - tan^^v^; 

Tgo3^ - Tcoi = tan^6'i^ sin2/3. 



(A.39a) 
(A.39b) 



Inserting this into eqs. (|A.38|) gives: 



' 2 Ml fi 



X 2 ,J Ml ' 



(A.40a) 
(A.40b) 



While perhaps not immediately obvious, it can be shown that eqs. (|A.40|) are indeed 
identical to eqs. ( [7071) 0. 

The fact that the Goldstone amplitudes (|A.40|) automatically yield more compact 
expressions shows that they can not only be used to check some amplitudes of Appendix 
A, but they also allow to derive the high-energy limit of the amplitudes for the production 



*Up to an overall sign between the ^Pq amplitudes; however, as shown in the Appendix of ref. the 
signs of noninterfering amplitudes are ambiguous when Majorana spinors are involved. 
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of longitudinal gauge bosons more quickly. We saw already that the amplitudes for the 
production of transverse gauge bosons vanish with some power of Mz / in the limit of 
large \m^\ if X is bino-like; the longitudinal helicity states therefore dominate the total 
production of gauge bosons in this important special case. (This is not true for higgsino- 
like or mixed LSP, however.) Most amplitudes for the production of Goldstone bosons 
can be read off directly from our results for the production of pseudoscalar and charged 
Higgs bosons, replacing sin/3 — >• — cos/3, cos/3 — >• sin/5 in the corresponding couplings of 
Higgs bosons to neutralinos and charginos. The only exception is H^G~ production, 
which gives the high-energy limit of H^W£ production. We list the nonvanishing partial 
wave amplitudes for completeness: 



: 2V2gWQ',^ 



-TO 



^/6 . 
~2 



vsm2/3g\Q',f-Q',f 



1 



1 



8 2(l + i?f) ^3(l + i?f)2 



1 + Rf 3(1 + Rff 
Rt 



l + Rl 



l + Rl 



3d " -0/3 2^0j Voj 

: j=vsm2ljg — 



2 

73^ 



1 + Rff 



(A.41a) 



(A.41b) 
(A.41c) 

(A.41d) 



In eqs.( |A.4T ), we have again only kept terms that remain finite at high energies. 
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Figure Captions 



Fig. 1 The rescaled LSP relic density Qh"^ as a function of the LSP mass for scalar 
SUSY breaking mass parameter m = 300 GeV and four combinations of the ratio 
tan/? of Higgs vevs and the mass rrit of the top quark. The solid (dashed) curves 
are for positive (negative) gaugino mass parameter M. Notice that M, the trilinear 
soft breaking parameter A and the supersymmetric Higgs (ino) mass parameter fi 
all vary along the curves, due to the relations between model parameters implied 
by radiative gauge symmetry breaking; in particular, the locations of the pole and 
thresholds caused by the light scalar Higgs boson are different in each case. 

Fig. 2 An enlargement of the threshold region of the solid curves in figs. la,b. For rrit = 
130 GeV (solid curve) the hh and Zh final states are accessible for > 75 and 
84 GeV, respectively, while for rrit = 160 GeV (dashed), they are only accessible for 
niy^ > 95 and 93 GeV, respectively. 

Fig. 3 The rescaled LSP density as a function of the LSP mass for two examples of higgsino- 
like LSPs. M, m and rrit are fixed for each curve, while /i, A and tan/3 vary. For 
M = —1 TeV (dashed curve) the LSP becomes bino-like at ~ 430 GeV, which 
explains the rapid rise of this curve at large m^. Only xx annihilation has been 
included in the calculation; this substantially overestimates Qh"^ for < Mw, as 
argued in the text. 

Fig. 4 Demonstration of the importance of the SUGRA imposed relations between particle 
masses. The solid curve shows the full SUGRA prediction for Qh"^ as a function 
of tan/3; eg. (p5|) implies that = mp/2 at tan/3=35.3, where P denotes the 
pseudoscalar Higgs boson. The long dashed curve has been obtained by keeping mp 
constant at the value SUGRA predicts for tan/?=2, i.e. mp = 783 GeV, but the 
f and b masses are still allowed to vary with tan/3. For the short dashed curve all 
effects of the b and r Yukawa couplings have been switched off. 

Fig. 5 Contours of constant Qh"^ in the plane spanned by M and tan/5, for fixed m = 250 
GeV and rrit = 140 GeV; results for all 4 combinations of signs of M and tan/? 
are shown. The solid and long dashed lines are contours where Qh'^=l and 0.25, 
respectively. The short dashed lines in (a) are contours where Qh^=0.025; since 
these contours are very close to the Z, h or P poles, in (b) - (d) we have merely 
indicated the location of these poles by the short dashed lines. The region outside of 
the dotted curves is exlcuded by various experimental and theoretical constraints, 
as described in the text. A flat universe requires 0.25 < Qh"^ < 1, while LSPs can 
build up the dark matter halo of galaxies if flh'^ > 0.025. 
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Fig. 6 Contours of constant Qh"^ in the plane spanned by M and tan/?, for the quadrant 
M > 0, tan/3 > 0, and four different combinations of m and rrit as indicated. The 
notation is as in figs. 5 (b)-(d). 

Fig. 7 The upper bound on \M\ which follows from the requirement that LSPs do not 
overdose the universe, fth'^ < 1. The long dashed, solid, and short dashed curves 
are valid for rrit — 110, 140 and 170 GeV, respectively, and show the absolute upper 
bound over the entire experimentally and theoretically allowed parameter space. 
The dotted curves emerge when one restricts the parameter space to |tan/3| < 15; 
the upper (lower) dotted curve is for rrit = 140 (170) GeV. In all cases the maximal 
|M| is found in the quadrant M < 0, tan/3 < 0. 

Fig. 8 Upper bounds on flh^ as a function of the mass mj^ of the SU (2) singlet sleptons. 
The long dashed, solid, and short dashed curves are for rrit — HO, 140 and 170 GeV, 
respectively, and present experimental bounds on sparticle masses; the long-short 
dashed curve is valid if the lower bounds on the masses of charged sparticles and 
gluinos are raised to 80 and 200 GeV, respectively, for rrit = 140 GeV. The dotted 
curve is based on a simple analytical approximation that assumes that the LSP is 
a pure bino, as described in the text. 
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